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Abstract 



We study the quantum transport through entropic barriers induced 
by hardwall constrictions of hyperboloidal shape in two and three spa- 
tial dimensions. Using the separability of the Schrodinger equation and 
the classical equations of motion for these geometries we study in detail 
the quantum transmission probabilities and the associated quantum res- 
f*^ . onances, and relate them to the classical phase structures which govern 

the transport through the constrictions. These classical phase structures 
are compared to the analogous structures which, as has been shown only 
recently, govern reaction type dynamics in smooth systems. Although the 
systems studied in this paper are special due their separability they can 
be taken as a guide to study entropic barriers resulting from constriction 
' geometries that lead to non-separable dynamics. 
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^ ■ 1 Introduction 

^ ■ A system displays reaction type dynamics if its phase space possesses bottleneck 

type structures. Such a system spends a long time in one phase space region 
}J] ■ (the region of 'reactants') and occasionally finds its way through a bottleneck 

5^ \ to another phase space region (the region of 'products') or vice versa. This 

type of dynamics does not only characterize chemical reactions but is of great 
significance in many different fields of physics and biology. Examples include 
ballistic electron transport problems [IJ, surface migration of atoms in solid state 
physics |2j , ionisation of Rydberg atoms in electromagnetic fields [31 14] , and on a 
macroscopic scale, the capture of moons near giant planets and asteroid motion 

In systems where the dynamics is smooth and Hamiltonian, the phase space 
bottlenecks eluded to above are induced by saddle-centre-. . .-center type equi- 
librium points, i.e. equilibrium points at which the matrix associated with 
the linearization of Hamilton's equations has one pair of real eigenvalues, ±A, 
and otherwise purely imaginary eigenvalues ±iwfc, k — 2, . . . , /, where / is the 
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number of degrees of freedom. In chemistry terms there is a 'transition state' 
associated with the bottleneck, i.e. a state the system has to pass 'through' 
on its way from reactants to products. The most efficient and commonly used 
approach to compute reaction rates is transition state theory, where the main 
idea is to place a dividing surface in the transition state region and compute the 
reaction rate from the flux through the dividing surface (for recent references, 
see the perspective paper [7]). This approach has major computational benefits 
over other methods to compute the reaction rate because the latter typically 
require the integration of trajectories in order to decide whether they are reac- 
tive (i.e. extend from reactants to products, or vice versa) or nonreactive (i.e. 
stay in the regions of products or reactants). Rather than this global informa- 
tion about trajectories, which to obtain is computationally expensive, transition 
state theory requires only local information about the phase space structures 
near the saddle-centre-., .-center equilibrium point - namely the construction 
of the dividing surface. However, in order to be useful and not to overestimate 
the reaction rate the dividing surface needs to have the property that it divides 
the phase space into a reactants and a products region in such a way that it 
is crossed exactly once by reactive trajectories and not crossed at all by non- 
reactive trajectories. The question how to construct such a dividing surface 
for systems with an arbitrary number of degrees of freedom has posed a major 
problem for many years, and has been solved only recently based on ideas from 
dynamical systems theory (see [1] and the recent review paper [S] with the ref- 
erences therein). The main building block in this construction is formed by a 
so called normally hyperbolic invariant manifold (NHIM) which is a manifold 
that is invariant under the dynamics (i.e. trajectories with initial conditions 
in the manifold stay in the manifold for all time) and is unstable in the sense 
that the expansion and contraction rates associated with the directions tangent 
to the manifold are dominated by those expansion and contraction rates as- 
sociated with the directions transverse to the manifold [S]. The NHIM is the 
mathematical manifestation of the transition state. In fact, the NHIM which 
is a sphere of dimension 2/ — 3 (with / again denoting the number of degrees 
of freedom) can be viewed to form the equator of the dividing surface which 
itself is a sphere of dimension 2/ — 2 located in a (2/ — l)-dimensional energy 
surface if it has an energy slightly above the energy of the equilibrium point. 
The NHIM separates the dividing surface into two hemispheres. All forward 
reactive trajectories (trajectories evolving from reactants to products) cross one 
of these hemispheres; all backward reactive trajectories (trajectories evolving 
from products to reactants) intersect the other of these hemispheres. More- 
over, the NHIM has stable and unstable manifolds. These have the structure 
of spherical cylinders R x S^^~^. Since they are of one dimension less than 
the energy surface they have sufhcient dimensionality to serve as impenetrable 
barriers in phase space [3] . They enclose the regions in the energy surface which 
contain the reactive trajectories and this way form the phase space conduits for 
reactions. 

Due to the spatial confinement, quantum effects are particularly strong for 
the passage through a phase space bottleneck, and accordingly, there is a strong 
interest in the quantum mechanical manifestation of the transition state. In 
molecular collision experiments, for example, high resolution spectroscopic tech- 
niques have been developed to directly or indirectly probe the transition state 
(see, e.g., |10|). Two quantum mechanical imprints of the transition state are 
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given by the quantization of the so-called cumulative reaction probability which 
is the quantum analogue of the classical flux, and the quantum resonances as- 
sociated with the transition state. The quantization of the cumulative reaction 
probability concerns the stepwise increase of the cumulative reaction probabil- 
ity each time a new transition channel opens as energy is increased. While this 
is quite difficult to observe in chemical reactions (see, e.g., the controversial 
experiment on the isomerization of ketene [TT]) this effect can be seen almost 
routinely as a quantization of the conductance in the ballistic electron transmis- 
sion through point contacts in semiconductor hetero-structures [12|, I13j , metal 
nano-wires [HI [15] and even liquid metals. The quantum resonances on the 
other hand, describe how wavepackets initialised on the transition state decay 
in time (see j8j for a detailed study). 

Moreover, there is a strong interest in the development of a quantum version 
of transition state theory, i.e. in a method to compute quantum reaction rates 
in such a way that it has similar computational benefits as (classical) transition 
state theory. Though much effort has been devoted to this problem it is still 
considered an open problem in the recent perspective paper [^. One major 
problem here seemed to be the lacking geometric insight which ultimately led 
to the realization of classical transition state theory. In [TBI [5] a quantum 
version of transition state theory has been developed which incorporates the 
classical phase space structures mentioned above in a natural way. It has been 
demonstrated to yield quite efHcient procedures to compute cumulative reaction 
probabilities as well as resonances. 

In this paper we are concerned with phase space bottlenecks which are not 
induced by equilibrium points. In the chemistry literature such bottlenecks are 
referred to as entropic barriers: in the microcanonical pictures this means that 
despite of the absence of a potential barrier, there is a minimum in the number 
(or to be more precise phase space volume) of possible configurations transverse 
to a reaction path. More concretely, we will consider potentialless systems 
with two and three degrees of freedom where the entropic barriers result from 
hard wall constrictions with the shape of an hyperbola and an (asymmetric) 
hyperboloid, respectively. We will be particularly interested in the phase space 
structures which govern the reaction dynamics in these systems and thus play an 
analogous role as in the case of a smooth Hamiltonian system with reaction type 
dynamics as mentioned above, and their quantum mechanical manifestations. 

The motivation for studying hyperboloidal geometries is that the resulting 
classical and quantum mechanical dynamics in such geometries are separable 
and in this sense completely solvable for such systems. This leads, as we will 
see, to a very transparent study of the influence of the phase space structures 
on the quantum transmission, and this way can serve as a flrst guide to study 
also non-separable dynamics in other constriction geometries. 

This paper is organized as follows. In Sec. [D we introduce in detail the 
systems studied in this paper and the associated transmission problems. In 
Sec. [3|we show how the Schrodinger equation of the transmission problem can 
be separated. The corresponding separations of the classical equations of motion 
are studied in Sec. U) The quantum and classical transmission probabilities are 
computed in Sec. [H Finally we compute and discuss quantum resonances in 
Sec. [6l and give a summary of the results and an outlook in Sec. [Tj 
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2 The transmission problem 



In the 2D case we consider a point particle moving freely in a region of the plane 
defined by 

9 2 

where {x, y) are Cartesian coordinates in the plane, and a and b are positive 
constants. We assume that, classically, the particle is specularly reflected when 
it hits either of the branches of the boundary hyperbola 

(see Fig. [ij . Quantum mechanically, this leads to the boundary condition that 
the wavefunction which describes the position of the point particle has to vanish 
on the boundary hyperbola ^ . In the wide-narrow- wide geometry of the region 
([1]) we can associate the part which has a; <C — 1 with the region representing the 
'reactants' and the part which has a; ^ 1 as the 'products', and that a 'reaction' 
has taken place when the particle has moved from reactants to products. This 
interpretation directly applies to the ballistic transmission of electrons through 
a point contact formed by a lead of the shape ((T]), but more generally can be 
viewed as a model describing the collective motion of a many body problem like 
a molecule from one configuration (or 'isomer') to another. 

In the 3D case we consider an analogous region in the three-dimensional 
space defined by 

222 

where {x, y, z) are Cartesian coordinates, and a, b and c are positive constants for 
which we impose the condition b > c. Note that this condition is only imposed 
for convenience and does not restrict the generality since one can simply swap 
the y axis with the z axis. The region ([3]) is bounded by the (asymmetric) 
hyperboloid 

222 
X y z^ 

- - + f- + - = 1 4 
52 

(see Fig. [T]). We again assume that, classically, the particle is specularly re- 
flected when it hits the boundary hyperboloid and hence also that the quantum 
mechanical (position) wavefunction vanishes on the boundary hyperboloid ([3]). 
We note that the region ([T]) in 2D can be formally obtained from the region 
^ in 3D by letting c — > which implies z — + 0. While taking this limit leads 
to no problems for the classical dynamics, one has to be more careful, due the 
Heisenberg uncertainty relation, when considering this limit in the quantum 
case. One can view the 2D transmission problem to be contained in the 3D 
transmission problem either by considering a small but finite c > which leads 
to a flat region near the x — y plane where for the energies under consideration 
no excitations in the z direction are possible, or by considering a cylindrical 
region in 3D where the base of the cylinder has the shape ([T]). 

The region ([1]) has a "bottleneck" contained in the y axis which is given by 
the line segment with minimal and maximal y values —b and -1-6, respectively. 
Similarly, the region ([3]) has a bottleneck in the y — z plane which is bounded 
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by the ellipse / h"^ + / c? = 1. In order to reduce the number of (effective) 
parameters we use as the length scale the maximum value of y in the bottle- 
neck. So formally we have b = 1 and the number of parameters specifying the 
accessible regions is 1 in the 2D case and 2 in the 3D case. 

The transmission through the bottlenecks can be viewed as a scattering 
problem. To this end we assume that a beam of (noninteracting) particles is 
incident from x <C — 1 (the 'reactants') and we want to compute the transmis- 
sion probability to a; ^ 1 (the 'products'). We will compute the transmission 
probability both classically and quantum mechanically in the spirit of transition 
state theory in Sec. [5l 

As mentioned in the introduction the motivation for choosing constrictions 
of the types ^ and ([U is that they are the most general type of hard wall 
constrictions for which the transmission problem can be separated and in this 
sense solved explicitly. We will discuss the separation in the following section 
(Sec. 13]) • In fact, in the 2D case the transmission problem is still separable if 
the constriction is composed of two branches of different confocal hyperbolas. 
However, the asymmetric case has no 3D analogue and we therefore restrict 
ourselves to the symmetric case ^ . Some aspects of the quantum transmission 
and the associated resonances through constrictions of the types and ^ 
have been addressed already in earlier papers. The quantum resonances for 
an asymmetric 2D constriction consisting of the branches of different hyperbola 
have been studied by Whelan [T7] . The quantum transmission problem (without 
resonances) through a constriction of the type ^ has been studied by Yosefin 
and Kaveh [T5]. Similarly, the transmission problem (again without resonances) 
has been studied for an axially symmetric hyperboloidal constriction in 3D by 
Torres, Pascual and Saenz ^9\, and for the asymmetric case by Waalkens [20] . 
The main purpose of the present paper is to study the quantum transmission 
and the assoicated resonances through the 2D and 3D constrictions ^ and ^ 
in a coherent way using the perspective of transition state theory. 

3 Separation of the Schrodinger equation 

For the quantum transmission problem, we have to find solutions of the free 
Schrodinger or Helmholtz equations 

or 

^2 / ^2 Q2 Q2 \ 

which for a; ^ 1 are waves propagating in the positive x direction and fulfill 
Dirichlet boundary conditions, i.e. we require the restriction of V' on the bound- 
ary hyperbola ([2]) resp. hyperboloid ([4]) to vanish. The Helmholtz equations 
([5]) and ([6]) together with their boundary conditions can be separated in elliptic 
and ellipsoidal coordinates, respectively, as we will discuss in the following two 
subsections which separately consider the 2D case and 3D case. 
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Figure 1: Accessible region confined by the boundary hyperboloid The 
region has a "bottleneck" in the y — z plane with the shape of an ellipse with 
semimajor axis 1 and semiminor axis c = — (in scaled coordinates). For 
the 2D case {c — 0, or equivalently c = 1), the accessible region is the area 
between the two branches of the hyperbola Q in the x — y plane. 
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Figure 2: Coordinate lines ^ =const. (dashed ellipses) and C =const. (solid 
hyperbolae) with the boundary hyperbola ^ in bold. The bold dots mark the 
focus points {x,y) = (0,±a). [a? = 5.) 
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3.1 The 2D system 

The 2D Helmholtz equation ^ together with the Dirichlet boundary conditions 
can be separated in eUiptic coordinates (C, £,) [1I1I121- Each of them parametrizes 
a family of confocal quadrics 



^ = 1- (7) 



where s € {C, ^} and = 1 + c? jl? . 

For s = ^ > a, both terms on the left hand side of Eq. ([7|) are positive and 
the equation defines a family of confocal ellipses with foci at = (0,±a). 

Their intersections with the x axis and y axis are at x = ^\/^ — a? and y = 
respectively. For a > s = C, > 0, the first term on the left hand side of 
Eq. ((7]) is negative giving confocal (two sheeted) hyperbolae with foci also at 
{x,y) = (0,±a). Their intersections with the y axis are at y = ±C; they do not 
intersect the x axis. 

The coordinate lines of C and r\ are shown in Fig. [2] Inverting Eq. ([7]) within 
the positive x — y quadrant gives 



. . (8) 
a 

y = (9) 

a 

with 

< C < a < (10) 

The remaining quadrants are obtained from appropriate reflections. However, 
it is also useful to reduce the discrete reflection symmetry of the system about 
the X axis and the y axis. In fact the solutions of the Helmholtz equation ([S]) 
fulfilling the Dirichlet boundary conditions along the boundary hyperbola ^ 
can be classified in terms of their parities and tTj, which correspond to the 
reflections about the y axis and x axis, respectively. We therefore introduce 
the symmetry reduced system which only has the positive x — y quadrant as 
the fundamental domain and impose Dirichlet (negative parity) or Neumann 
boundary conditions (positive parity) on the x and y axes. The Cartesian 
coordinate axes are obtained from the elliptic coordinates C and ^ in terms of 
the equalities in (fTU|) : C = gives the x axis; — a gives the segment of the y 
axis between the focus points, the rest of the y axis has C = a (see Fig.[31[a)). 

The boundary hyperbola ([2|) (in scaled coordinates) coincides with the co- 
ordinate line C = li i-e. in the region ([T|) Q takes values in [0, 1]. Considering 
only the region enclosed by the boundary hyperbola ([2]), the coordinate lines 
^ =const. > a are transverse to the x direction. To this end note that the sin- 
gular coordinate line ^ = a contains the 'bottleneck' (x, y) e {0} x [—1, 1]. The 
coordinate ^ thus parametrizes the direction of the transmission; ^ parametrizes 
the direction transverse to the transmission. 

The parameter a determines how strong the narrowness of the constriction 
changes with x: for a — > oo the constriction becomes an inflnitcly long rcctan- 
gualar strip; for a — )■ 1 the constriction degenerates to the y axis with a hole of 
width 2 about the origin. 
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Figure 3: (a) Singular elliptic coordinate surfaces, for the 2D system. The line 
segment [0, 1] is half of the "bottleneck" on the y axis, (b) Singular ellipsoidal 
coordinate surfaces, for the 3D system, (c) "Bottleneck" (shaded region) in the 
y — z plane bounded by the 3D hyperboloidal constriction whose intersection 
with the y — z plane is the ellipse + z'^/{l — c^) — 1, and singular coordinate 
patches ^ = a and r] = a (inside and outside of the eUipse y^/a'^ + z^/ (a^ — c^) = 
1, respectively). 

With the ansatz '0(CjC) = '0c(C)'0c(O the partial differential equation ([5]) 
can be separated and turned into the set of ordinary differential equations 

-^(^^/^^^Jms)^e{s'~sI)Ms), (11) 

where s G {(,£.} and s| denotes the separation constant. The equations for ( 
and ^ are identical, but they have to be considered on the different intervals 
(jlOp and for different boundary conditions. In fact the equations have regular 
singular points |21| at ±a. These regular singular points have indices and 1/2, 
i.e. there are solutions, which near ±a are of the form iJjs{s) = {s T a)'^ip{s) 
where V'(s) is analytic and q — or q = 1/2. As the elliptic coordinates 
give for the regular singular point a the Cartesian y axis, the indices determine 
the parities tt^ of the total wave function ip{C,£,) [23; i-^- q = or q = 1/2 
correspond to total wave functions which have — + or tt^ — respectively. 
The value of "0^ at the ordinary point C = determines the parity iTy . 

For the computation and interpretation of the results below, it is useful to 
remove the singularities in (|lip . This can be achieved by the transformation 

(C(^),e(A)) =a(cos(i^),cosh(A)), (12) 

which is the standard parametrization of elliptic coordinates by triginometric 
functions. Inserting (fT^ into ^ and ^ gives 

X = a sin(i^) sinh(A) , y = a cos(z/) cosh(A) . (13) 

To cover the positive x — y quadrant {v, A) have to vary in the intervals 

0<i/<7r/2, 0<A<oo. (14) 

The boundary hyperbola dH) has 

v = vb ~ arccos(l/a) . (15) 

Extending the intervals (fH]) to 

vb l£ l£ ''^ — i^B , —oo < X < oo . (16) 
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we get a full regular cover of the region ([T]) in terms of the strip [i^b, tt — j^b] x H. 
Transforming to the coordinates {ly, A) leads to 

- ^^M^) = '^^E {s'is) - si) Ms) , (17) 

where s S {i^, A}, s(s) S {CC*^); ^l-^)} a-rc the functions from (112p and the Us are 
the signs cr\ = + and ctj/ = — . Each of these equations can be interpreted as a 
one-dimensional Schrodinger equation with a Hamiltonian of the standard type 
H = —h?{d'^ /dx'^)/2 + V ("kinetic plus potential energy") with effective energy 
and potential 

Es,eS = -asEsj , VsMs) = -(TsEs^{s) . (18) 

The effective energies and potentials (fTB]) are shown for "representative" values 
of the separation constant in Fig. [Tl^a) of Sec. ID Here v varies in an interval 
of length TT which is the period of the effective potential Vj^^off • What we mean by 
"representative" will be explained in Sec.[31 where we analyze the corresponding 
classical system. Since the effective potential V\,es is symmetric under the 
reflection A i— > —A (i.e. the reflection about A = 0) there are solutions of 
P7p that are symmetric or antisymmetric under this reflection. Using (|13|) 
we can relate the behaviour of solutions under this reflection to the parity tt^^. 
Similarly, since the effective potential V^^cS is symmetric under the reflection 
1/ 1-^ —v + TT (i.e. the reflection about i/ — 7r/2) there are solutions of (fT7|) 
that are symmetric or antisymmetric under this reflection. Again using (|13p 
we can relate the behaviour of solutions under this reflection to the parity iTy. 
The parities tt^ and TTy are marked at the top of Fig. [7Ja). The fact that the 
algebraic counterparts of ([T7]) in ([TT|) are identical, is reflected in ([T7| by the 
substitution v ^ iv which relates the equation for v to the equation for A. 



3.2 The 3D system 

Similarly to the 2D case, the 3D Helmholtz equation ^ together with the 
Dirichlet boundary conditions can be separated in ellipsoidal coordinates (C, ry, ^) 
[2T| [22] . Each of them parametrizes a family of confocal quadrics 

2 9 9 

y 

72 12 72 l2 ^ ' 



where s € {CiVtOj = 1 — 2^/6^ and = 1 + a^/&^. 

For s = ^ > a, all terms on the left hand side of Eq. (I19p are positive and 
the equation deflnes a family of confocal ellipsoids. Their intersections with the 
y — z plane, the x — y plane and the x — z plane are planar ellipses with foci at 
(y, z) — (±c, 0), {x, y) = (0, ±a) and (x, z) = (0, ±(a^ — c^)-^/^), respectively. For 
a > s — rj > the flrst term on the left hand side of Eq. (fT^ becomes negative. 
Eq. (fT9|) thus gives confocal one sheeted hyperboloids. Their intersections with 
the y — z plane are planar ellipses with foci at (y, z) — (±c, 0); the intersections 
with the X — y plane and the x ~ z plane are planar hyperbolas with foci at 
(x, y) — (0, ±a) and (x, z) = (0, ±(a^ — c^)^/^), respectively. For c> s = C, > Q, 
the flrst and third terms on the left hand side of Eq. IT^ are negative giving 
confocal two sheeted hyperboloids. Their intersections with the y ~ z plane 
and the x — y plane are planar hyperbolas with foci at [y, z) — (±c, 0) and 
(x, y) — (0, ±a), respectively; they do not intersect the x — z plane. 
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Figure 4: Ellipsoidal coordinates surfaces ^ = 5 (ellipsoid in red), ry = 1 (the one 
sheeted boundary hyperboloid (H]) in blue) and C, = 1/2 (two sheeted hyperboloid 
in green), for (a^,c^) = (5,0.2). 



The coordinate surfaces of C, r] and ^ are shown in Fig. H) Inverting Eq. (fTO|) 
within the positive x — y — z octant gives 



Vie 


-a2)(a2 _ 










eK 






ac ' 










C2)(c2-C') 




cV a? — 


C2 



(20) 
(21) 
(22) 



with 



0<C<c<?/<a<^. (23) 



The remaining octants are obtained from appropriate reflections. Again, we 
also introduce a symmetry reduced system which has the positive x — y — z 
octant as the fundamental domain. The solutions of the Helmholtz equation © 
fulfilling Dirichlet boundary conditions along the boundary hyperboloid ^ with 
parities tt^, tTj, and tt^ are then obtained from the symmetry reduced system 
by imposing Dirichlet or Neumann boundary conditions along the Cartesian 
coordinate planes which in terms of the elliptic coordinates (C, r], A) are given 
by the equalities in one of the equations in (^5]) : C = gives the x — z plane; 
C — c and rj = c give two surface patches which together cover the x ^ y plane; 
•q = a and £, = a give two surface patches which together cover the y — z plane 
(see Fig.[3i;b)). 

The boundary hyperboloid @ (in scaled coordinates) coincides with the 
coordinate surface 77 = 1, i.e. within the region ([3]) is restricted to [c, 1]. 
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Considering only the region enclosed by the boundary hyperboloid the co- 
ordinate planes ^ =const.> a are transverse to the x direction. Note that the 
singular coordinate plane ^ = a is a region in the y — z plane which is en- 
closed by an ellipse which lies outside of the hyperboloidal constriction (see 
Fig. Eljc)). The coordinate <^ thus parametrizes the direction of transmission; rj 
and ^ parametrize the two directions transverse to transmission. 

The parameter c determines the asymmetry of the cross-section of the con- 
striction with c = leading to an axially symmetric constriction and c = 1 
leading to the 2D case. The parameter a determines how strong the narrowness 
changes with x: for a — > oo the constriction becomes cylindrical with an ellipti- 
cal cross-section; for a 1 the constriction degenerates to the y — z plane with 
a hole having the shape of an ellipse. 

With the ansatz ipiCrj,^) — '0c(C)V'r)('7)'/'5(^) the Helmholtz equation ® 
can be separated and turned into the set of ordinary differential equations 

- 1^ (VF^^^K^^^^) Ms)^E{s^-2ks' + l)Ms), (24) 

where s G {C,Vt^} ^'^d ^ ^-nd I denote the separation constants. The equations 
for ^, 1] and ^ are identical, but they have to be considered on the different 
intervals ([23|l and for different boundary conditions. For later purposes it is 
useful to rewrite (1241) in the form 



(^^(s2-a2)(s2_c2)A^ ij,is)^E{s'-sl){s'-sl)Ms), (25) 

sl^k + {e-lf'\ (26) 



2m 

where 



si 



and conversely 



k^\{sl + sl), l = slsl. (27) 

Similarly to equations (fTTj) in the 2D case the equations (p4|) have regular 
singular points [21] at ±a and ±c. All these regular singular points again have 
indices and 1/2 like in the 2D case. Thus there are solutions, which near 
a — ±a or a = ±c are of the form ipsis) = (s — a)'^'"ip{s) where ip{s) is analytic 
and Qa — or qa = 1/2. As the ellipsoidal coordinates {CVtO gi'^^ foi" the 
regular singular points ±a and ±c the Cartesian y — z plane and x — y plane, 
respectively, the indices determine the parities tTx and tt^ of the total wave 
function i'{C,V:0 E2- More precisely, qa = or qa — 1/2 correspond to total 
wave functions which have tt^ = + or — —, respectively, and qc = or 
qc = 1/2 correspond to total wave functions which have tt^ = -I- or tt^ = — , 
respectively. As in the 2D case, the value of 4'c ^t the ordinary point C = 
determines the parity ny. 

For the computation and interpretation of the results below it is useful to 
remove the singularities in This can be achieved by the transformation 



{({,,), r7(/i), e(A)) - a{q sn(^., q), dn{pi, q'), ^^^) , (28) 

cn(A,g) 

where sn(0, (7), cn(0, g) and dn{(f>,q) are Jacobi's elliptic functions with "angle" 
(f) and modulus q [23] . Here the modulus is given hy q — c/a and g' = (1 — 
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gr2-)i/2 (igjiotes the conjugate modulus. This is the standard parametrization of 
elhpsoidal coordinates by eUiptic functions [U]. 

Expressing the Cartesian coordinates in terms of {iy,fi, A) gives 

, sn{X,q)sn{fi,q')dn{iy,q) 

x = q a — ^ , 

cn(A, q) 

dn(A,(7)dn(^,g')sn(t/,g) 

y = a 7T — ^ . (29) 

cn(A,q) 

, cn(/i, q') cn{iy, q) 

z = q a — r . 

cn(A, q) 

To cover the positive x — y — z octant (i^, /i, A) have to vary in the intervals 

0<v< K{q) , 0<fi< K{q') , < A < K{q) , (30) 

where K{q) and K{q') are Legendre's complete elliptic integral of first kind with 
modulus q and q', respectively. The boundary hyperboloid ([IJ has 

M = /iB = F{{{a^ - l)/(a2 - c^)f'\ q') , (31) 



where F is Legendre's incomplete elliptic integral of first kind which in ([3T|l has 
argument ((a^ — l)/(a^ — c^))^/^ and modulus q' . Extending the intervals ([501) 
to 

0<iy<4K{q), MB < M < 2if(9') - MB , -Kiq)<\<Kiq) (32) 

we get a double cover of the region ^ in terms of the 'soHd torus' R/ {4:K{q)1i) x 
[HB,'2K{q') - Mb] x [-K{q), K{q)], where R/{AK{q)Z) denotes the topological 
circle resulting from identifying points in R differing by integer multiples of the 
period in ly which is 4K(q). In Fig. [5] we present the solid torus as the cube 
((5^ . where the opposite sides j/ = and i/ = 4:K{q) have to be identified. Each 
of the smaller cubes 

[n,K{q), {n, + l)Kiq)] x [0, K{q')] x [0, ±K{q)] and 
KA'(q), {n, + l)K{q)] x [K{q'),2K{q')] x [0, ±Kiq)] 

in Fig. El with Cz "Zi represents one Cartesian x — y — z octant of the region pl). 
Note that each of the smaller cubes ([551) has four neighbours. This property 
can be understood from the fact that in order to regularise the coordinates 
(C,?7,^) in terms of the coordinates (z^, /x, A) we have to regularise each of the 
four singular transition between two x — y—z octants shown in Fig. [31 (note 
that the singular patch 77 = a in Fig. [3lis not accessible in The two covers 
of the double cover are related by the involution 

S{y, M, A) - {2K{q) 2K{q'),X) , (34) 

which leaves the Cartesian coordinates (12^ fixed (see also Fig. (51). 
Transforming (|24p to the coordinates (i^, /i, A) leads to 

- I^T^V-Ks) = - 2fcs2(s) + I) Us) , (35) 

2m ds^ 
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Figure 5: Representation of tlie solid torus ([32|) as a cube with periodic boundary 
conditions in v. Each sniaU cube represents one Cartesian x — y — z octant. The 
octants corresponding to the smaUer cubes are indicated by a 'binary' labeling 
with respect to the signs of x, y and z (e.g., (-,-,-) corresponds to 0, (-,-,+) 
corresponds to 1, etc.). The shaded planes mark the boundary hyperboloid 
which on the double cover is given by /i = /j,b and /i = 2K'{q) — /^b, where 
K'{q) = K{q'). 
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where s € {z^, /i, A}, s(s) £ {C{i'),r]{^),^{X)} are the functions from (^5]) and 
the (Tg are the signs a\ — = + and cr^ = — . Each of these equations can 
be interpreted as a one-dimensional Schrodinger equation with a Hamiltonian 
of the standard type H = — /i^(d^/da:^)/2 + V ("kinetic plus potential energy") 
with effective energy and potential 

Es.cS = 'Js^l, Vs.csiS) ^ ~<Js^{sHS) ~ 2ks\s)) . (36) 

The effective energies and potentials p6p are shown for representative (again see 
Sec. lU values of the separation constants k and I in Fig. fTOTa) of Sec. HI where 
/i and 1/ vary in intervals of length 2K{q') and 2K{q), which are the periods of 
the effective potentials V^i^eff and V^^eS, respectively. 

The reflection symmetry of the effective potential V\,cs about A = leads 
to solutions of psp that are symmetric or antisymmetric under this reflection. 
Similar to the 2D case we can use (|29|) to relate the behaviour of solutions 
under this reflection to the parity tt^. The effective potential V^,eff is symmetric 
about /i = K{q'), and using psp the symmetry or antisymmetry of solutions of 
((55)) under the corresponding reflection jj, i-^ 2K{q') — ^ can be related to the 
parity tt^. The effective potential V^^es has reflection symmetry about v ~ 
and V = K{q). Eq. ([29]) relates the symmetry or antisymmetry of the solutions 
under the corresponding reflections v i—^ —v and v i— > 2K{q) — v to parities 
TTy and TTz, respectively. We note that, like their algebraic counterparts (IM)) . 
the wave equations psp for v, fi and A are identical, if one considers them on 
different intervals (in the complex plane). The equations for fi and A can, e.g., 
be related to the equation for u using the identities sn{u + K{q) + iK{q'),q) = 
g""^ dn(u, g)/ cn(u, g) and sii{—m + K{q)+iK{q'),q) = q~^dn{u,q') in ([28]) . 
This is similar to the statement on the wave equations PT|) in the 2D case. 

4 The classical systems 

We will now study the classical dynamics of the transmission problem described 
in Sec. [21 As mentioned in Sec. [2] the classical motions consist of motions along 
straight lines in the regions ([T]) and ([3]) with specular reflections at the bound- 
ary hyperbola and hyperboloid, respectively. Like the Helmholtz equations with 
the Dirichlet boundary conditions imposed along the boundary hyperbola and 
hyperboloid the classical equations of motion can also be separated in elliptic 
(2D) and ellipsoidal coordinates (3D). The separability implies that the classical 
dynamics is integrable, i.e. there are as many constants of the motion (the sepa- 
ration constants) that are independent and in involution as degrees of freedom. 
A modiflcation of the Liouville- Arnold theorem (25j says that the space of the 
classical motion is (up to singular sets of measure zero) foliated by invariant 
cylinders (the analogues of invariant tori in closed systems). In the following we 
will have a closer look at these foliations for both the 2D and 3D system. 

4.1 The 2D system 
4.1.1 Phase space foliation 

Separating the equations of motions for the free motion in the plane in the 
elliptic coordinates {(,£,) introduced in Sec. 13.11 yields that the momenta ps 
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Figure 6: Bifurcation diagram in terms of the variable for the 2D system, 
with =3/2 (a), and in terms of the variables {si, S2) (b) and (fc, I) (c) for the 
3D system, with (o^c^) = (3/2,1/2). 

conjugate to s, s G {C^C}: are given by 

pI = 2mE^^ (37) 

(see f5^), where is a separation constant which acts as the square of the 
turning point of the respective degree of freedom s G {CjC}- These equations 
are the analogues of the separated Helmholtz equations in the algebraic form 
pip . Similarly, for the coordinates s € (j^, A) and their conjugate momenta p^, 
the analogue of the regularized separated Helmholtz equations (|17p are given by 

pj = as2mE{s^{s) - sj) (38) 
= 2m{Es,,s-VsMs)), (39) 

where s e {i', A} and s{s) € {C('^)jC(A)} in are the functions defined in 
(fT^ , and the effective energy and potential in ([S^ are defined as in (fTH|) . 

The specular reflection at the 2D boundary hyperbola ^ = 1 or equivalently 
1/ = i/B and 1/ = TT — i/B is described by mapping the phase space coordinates 
right before the reflection to the phase space coordinates right after the reflection 
according to 

or 

{iy,X,Pi,,Px)i-^{i^,X,~Pi,,Px), (41) 

respectively. As opposed to the phase space coordinates (s,Ps), s £ {CiS,} the 
phase space coordinates {s,ps), s G {f^, A} lead to a smooth description of the 
motion (apart from the specular reflections). 

Like in the quantum case in Sec. 13.11 we can also introduce a symmetry 
reduced system in the classical case. For the symmetry reduced system the 
motion is confined to the positive x — y quadrant of the region |T]) with specular 
reflections not only at the boundary hyperbola @ but also at the Cartesian 
coordinate axes. 
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Figure 7: (a) Effective potentials and energies for the types of motion Ti, T2 and 
T3 defined in Fig. Elja). For the v degree of freedom the hatched regions mark 
the forbidden regions [0, v-q\ and \k — v^^ tt] which are not contained in the region 
([T]). (b) Phase curves parametrized by s\ (with E — const. > 0). For s\ = 0, we 
have the forward and backward reaction paths (the free flight motions along the 
X axis) which correspond to the branches of the green solid hyperbola in the left 
panel and the green central dot in the right panel. For a fixed S2 ^ (0; 1)? we 
have two invariant cylinders of forward and backward reactive trajectories which 
do not involve specular reflections at the boundary hyperbola ^ . These appear 
as the two branches of the black solid 'horizontal' hyperbola in the left panel 
and the inner black solid circle in the right panel. The two cylinders which have 
S2 = 1 and have forward and backward trajectories that touch the boundary 
hyperbola ^ tangentially are marked by the red dashed curves in either panel. 
A fixed S2 G (Ij a^) represents two cylinders which have forward and backward 
reactive trajectories that involve specular reflections at the boundary hyperbola. 
These are marked by the two branches of the black solid 'vertical' hyperbola in 
the left panel and the corresponding chopped circle in the right panel. The value 
s| = represents the periodic orbit (or 'transition state') TS which corresponds 
to the origin in the left panel and the corresponding blue chopped circles in 
the right panel, and its stable and unstable manifolds with their forward and 
backward branches W^^^^ forming the blue cross in the left panel and coinciding 
with the blue chopped circles in the right panel. A fixed > represents two 
cylinders of nonreactive trajectories on the reactants side (A < 0) and products 
side (A > 0), respectively. These correspond to the two branches of the black 
solid 'vertical' hyperbola in the left panel and the corresponding chopped circle 
in the right panel. The black dashed lines in the left panel mark the forward 
{p\ > 0) and backward {p\ < 0) dividing surfaces DSb/f. In the right panel 
these appear as the hatched chopped disk, (a^ = 5.) 
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The physical meaning of the separation constant becomes more clear 
from multiplying it with (2mi?)^/^ and expressing it in terms of Cartesian co- 
ordinates. A little bit of algebra then gives 



V^sl = Ll + aVy = - {Ll_ + Ll+) , (42) 

where — xpy — ypx is the angular momentum about the origin, and L^^ — 
{x + a)py — ypx and Lz+ = {x — a)py — ypx are the angular momenta about 
the focus points {x,y) = (0,±a). This is the second constant of the motion 
beside the energy which makes the system integrable. A modification of the 
Liouville- Arnold theorem then implies that the four dimensional phase space is 
foliated by invariant cylinders (see below) which are given by the common level 
sets of the constants of motion E and or equivalently E and s\. In fact 
the energy plays no major role for the motions. It just determines the speed 
of the motion along the straight lines (in configuration space). As indicated by 
the occurrence of the energy as a multiplicative factor in the equations for the 
separated momenta ([BT)) and ([38]) . energy surfaces of different positive energies 
only differ by the scaling of the momenta, and accordingly they all have the 
same type of foliation by invariant cylinders. To discuss the foliations of the 
energy surfaces it is thus sufficient to consider a single energy surface of fixed 
energy > 0. The different types of cylinders contained in the energy surface 
of this energy are then parametrized by the second constant of motion, s|, in 
the following way. 

First of all, in order to simultaneously have real momenta in the physical 
ranges C S [0, a] and ^ € [a,cx3) (see (flU)) ') the separation constant can only 
take nonnegative values. We therefore will occasionally write S2 — \fs\- The 
interval [0, cxd) contains three subintervals which correspond to different smooth 
families of cylinders which we denote by 

Ti : < < 1 , 

T2 : \<s\<a^ , (43) 
T3 : <s\. 

At the values s\ = s\ — \ and s\ = a} the families of cylinders bifurcate, 
and these parameter values thus present critical motions to which we will come 
back below (also see the bifurcation diagram in Fig.[S]). 

To understand the motions on the different types of cylinders Ti, T2 and T3 
it is useful to consider the corresponding effective potentials and energies (fT8|) 
and phase portraits in the planes v — py and \ ^ p\ in Fig. [7] in combination 
with the projections of the cylinders to configuration space which are shown in 
Fig.H 

The common level set of the constants of motion E and s\ in Ti consists 
of two disjoint cylinders which both extend over all values of x. On one of 
these cylinders p\ is always greater than zero, and on the other p\ is always less 
than zero. These cylinders are thus foliated by forward and backward reactive 
trajectories, respectively. The motion oscillates in the v degree of freedom in 
such a way that the trajectories do not hit the boundary hyperbola The 
topology of the cylinders, E, x SS^, becomes apparent from taking the Cartesian 
product of the lines (~ R) in the phase plane \ — p\ in the left panel of Fig. [71(b) 
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with the corresponding topological circle in the phase plane v — in the right 
panel of Fig.[7i;b). 

Similarly, a common level set of the constant of motion in T2 consists of two 
disjoint cylinders of which one again consists of forward reactive trajectories and 
the other again consists of backward reactive, but the oscillations in the trans- 
verse degree of freedom u now involve reflections at the boundary hyperbola. 
To simplify the discussion, we will glue together the two line segments in the 
V — Pv plane in the right panel of Fig. [TJ^b) which have positive and negative p^, 
respectively, at the points v = and v — t: — i.e. ai v = and = tt — z^b 
we identify Pi, and —p^. Note that strictly speaking the momenta are not de- 
fined along the boundary hyperbola. However, the gluing can also be justified 
from physical considerations by viewing the hard wall potential which causes 
the reflections as the limiting case of a smooth potential that becomes steeper 
and steeper. The resulting object can then again be viewed as a topological 
circle, S"'^, and taking the Cartesian products with the corresponding lines ~ R 
in the X — p\ planes we again obtain topological cylinders R x similar to 
those in Ti. 

In contrast to the cylinders above, the common level set of the constants of 
motion in T3 consists of two disjoint cylinders which when projected to config- 
uration space are both bounded away from the y axis by the ellipse = s|. 
These cylinders are foliated by nonreactive trajectories which stay on the side 
of reactants and products, respectively. 

The critical value 53 = 1 corresponds to the limiting motion between Ti 
and T2. The level set of the constants of motion E > and S2 = 1 consists of 
two disjoint cylinders which contain forward and backward reactive trajectories, 
respectively, which hit the boundary hyperbola ^ tangentially (see the dotted 
line in the right panel of Fig.jTjb)). 

At the critical value S2 = the two cylinders of type Ti degenerate to two 
lines ^ R given by the Cartesian products of the dot at the centre of the v — p^ 
plane in the right panel of Fig. [T^b) with the corresponding lines in the X — p\ 
plane in the left panel of the same figure. One of these lines corresponds to a 
trajectory along the x axis which has px = (2mi?)^/^; the other line corresponds 
to a trajectory along the x axis which has Px — —{2mEy/'^. These can be 
viewed as the forward and backward reaction paths, i.e. they are the unique 
trajectories, which for a fixed energy, are reactive and do not involve any motion 
in the transverse degree of freedom (2^1 [5] . 

The critical value s| — c? represents the unstable periodic orbit along the y 
axis which, for a fixed energy i? > 0, bounces back and force between the two 
branches of the boundary hyperbola The common level set of > and 
s\ = a} consists not only of this periodic orbit but also of the stable and unstable 
manifolds W and of this periodic orbit. In the X~p\ phase plane in the left 
panel of Fig. [71(b) the stable and unstable manifolds occur as the cross shaped 
structure which has the periodic orbit at the center. With our interpretation 
of reflections at the boundary hyperbola to be smooth the periodic orbit has 
the topology S^, and its stable and unstable manifolds are cylinders R x S^. 
The stable and unstable manifolds are of special significance for the classical 
transmission since they are of codimension 1 in the energy surface, i.e. they 
have one dimension less than the energy surface, and this way have sufficient 
dimensionality to act as impenetrable barriers in the energy surface [27] . In fact, 
the stable and unstable manifolds form the separatrices between reactive and 
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nonreactive trajectories. More precisely, and W" each have two branches: 
we denote the branch of W which has p\ > (resp. p\ < 0) the forward (resp. 
backward) branch, Wj (resp. W^), of the stable manifold. Similarly we denote 
the branch of which has px > (resp. px < 0) the forward (resp. backward) 
branch, Wj^ (resp. W^), of the unstable manifold (see Fig. [Tl^b)). Moreover, we 
call the union of the forward banches, 

Wf ■.= W^UWf , (44) 

the forward reactive cylinder, and the union of the backward branches, 

Wb := W,' U , (45) 

the backward reactive cylinder. The forward reactive cylinder encloses all tra- 
jectories in an energy surface of the respective energy E > which are for- 
ward reactive; the backward reactive cylinder encloses all trajectories in such 
an energy surface which are backward reactive. The nonreactive trajectories 
are contained in the complement of these regions. The forward and backward 
reactive cylinders thus play a crucial role for the classification of trajectories 
with respect to their reactivity. They can be viewed to form the phase space 
conduits for reaction. In particular, the forward and backward reaction paths 
mentioned above can be viewed to form the centerlines of the regions enclosed 
by these cylinders. 

The periodic orbit, or more precisely, the family of periodic orbits oscillat- 
ing along the y axis with different energies E > can be viewed to form the 
transition state or activated complex. Reactive trajectories of a given energy 
E > pass 'through' the periodic orbit at that energy (the 'transition state at 
energy i?') in the following sense. Setting a; = on the energy surface defines a 
two-dimensional surface in the energy surface which is given by 

DS = {{x,y,Px,Py) e E,"* : x = 0, y e [-1,1], 
pI+pI = 2mE} 

4 (46) 
= {{v,\,Pty,px) e ]R : A = 0, e Wb,t^ - i^b], 

Px+pI^ '2mEa'^{l ~ cos^ v)} . 

With our convention to identify the momenta — and -l-pu ai v — vb and 
I' — TT — i^B, the surface DS has the topology of a two-dimensional sphere, S^. 
It defines a so called dividing surface that has all the desired properties that 
are crucial for the transition state computation of the classical transmission 
probability from the flux through a dividing surface. First of all, it divides the 
energy surface into a reactants part {x < 0) and a products part {x > 0). In 
order to be reactive a trajectory thus has to intersect the dividing surface. In 
fact the periodic orbit or transition state at energy E given by 

TS = {(x,2/,p„Py)eR2xR2 : x ^ 0, p^ ^ 0, y e [-1,1], pi = 2mE} (47) 



can be viewed to form the equator of the dividing surface (I46p . It separates the 
dividing surface into two hemispheres which we call the the forward dividing 
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Figure 8: Configuration space projections (shaded regions) of the invariant 
cylinders corresponding to the three types of 2D motions Ti, T2 and T3. The 
bold lines mark the boundary hyperbola The ranges for x and y are both 
[-3/2,3/2]. (a2 = 3/2.) 



surface 

DSf = {(a;,y,p^,Py) € : x = 0, y G [-1,1], 

pl+pl = 2mE,p,>0} ^^^^ 
= {{v,X,p„,px) e IR"* : A = 0, e [i^B.TT- j^b], 
pI+pI ^ 2mEa^{l - cos^ ly), px> 0} . 
and the backward dividing surface 

DSb = {{x,y,Px,Py) e R'^ : X = 0, y e [-1,1], 

pl+pl = 2mE,p,<0} ^^^^ 
= {{v,X,p^,px) e M"* : A = 0, 1/ e [^'B,7r- i'b], 
pI+pI = 2mEa^{l - cos^ i^), px < 0} ■ 

These two hemispheres appear in the right panel of Fig. [Tl^b) as the disk en- 
closed by the blue curve that represents the transition state periodic orbit TS 
in the v — Pi, plane. Note that the circles contained in this disk have to be com- 
bined with the two corresponding lines in the \ — px plane in the right panel of 
Fig. Efb) which have px > Q (corresponding to forward reactive trajectories) or 
Pa < (corresponding to backward reactive trajectories). All forward reactive 
trajectories have a single intersection with the forward dividing surface, and 
all backward reactive trajectories have a single intersection with the backward 
dividing surface. Nonreactive trajectories do not intersect the dividing surface 
at all. The dividing surface is everywhere transverse to the Hamiltonian flow 
apart from its equator, which is a periodic orbit and thus is invariant under the 
Hamiltonian flow. 



4.1.2 Action integrals 

In the previous subsection we have seen that the phase space is (up to criti- 
cal motions which form a set of measure zero) foliated by invariant cylinders 
where the cylinders are given by the Cartesian products of circles in [v^p^) and 
unbounded lines in {\,px). For the {v^Pu) component of the motions we can 
directly introduce action-angle variables As we will see below we can also 
associate an action type integral with the unbounded {\,px) component of the 
motion. Both of these actions will play a role in the semiclassical computation of 
the cumulative reaction probability and the quantum resonances (see Sections [S] 
andini respectively). 
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Action integrals depend on the type of motion, and typically change from one 
type of motion to another. For the actions associated with the v or equivalently 
C degree of freedom we find 

/c = i /pcdC = \l ^l^dc , (50) 

where we took from p7p and the integration boundaries C- and are given 
by C- =0 and = S2 for motions of type Ti and = 1 for motions of type T2 
and T3. The corresponding action integral for the symmetry reduced system, 
which we denote by is given by 

^C-^^C- (51) 

To understand the analytic nature of the action integral we substitute 
z = in Eq. ^ which gives 

Ir = V2^- r ^^^4^clz , (52) 

TT W{Z) 
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where 

w^{z)=P:,{z):=\{{z-z,), (53) 

i=l 

and z_, are consecutive elements of the set {zi = 0, Z2 = a^, Z3 = S2, Zb — 1}. 
Here Zb = 1 corresponds to the boundary hyperbola. The differential dz/w{z) 
has the four critical points {zi, Z2, z^, 00} which means that the integral (|52ll is 
elliptic. We refrain from expressing this integral in terms of Legendre's standard 
integrals [25] • Instead, and for later purposes (see Sections [5] and 0, we interpret 
the integral for motions of type T2 and T3 as an Abelian integral on the 
elliptic curve 

T^,^{{s,w)e€^ : w^^Psiz)}. (54) 

Here C defines the compactified complex plane (i.e. the Riemann sphere). The 
algebraic curve is of genus 1, i.e. it has the topology of a 1-torus. For 
motions of type T2 or T3 the action in (j53p can then be written as 

I^^V2^:^ f (z-4)-, (55) 

where for T2 , the integration path 7^ is defined as illustrated in Fig. [11(b) . For 
T3, the order of s| and along the real axis in Fig. [H^b) is reversed. However 
this does not affect the definition of 'y^ for T3. Due to the billiard boundary 
the integration path 7^ is not closed on F^,, i.e. the integral is an incomplete 
elliptic integral. 

On Fu, we also define the complete elliptic integral 

h^iV2^—f (z~sj)— (56) 
27r ./^, w 
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and its symmetry reduced partner 



k = hi (57) 
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Figure 9: (a) The graph of the polynomial P3 defined in ((SH) for a real separation 
constant s\ satisfying \ < s\ < (motion type T2). (b) Elliptic curve for 
= Pz{z) with P3 as in (a), and integration paths 7^ and 7^ which define the 
integrals in (|55|) and ([56| . Only one half (Riemann sheet) of the elliptic curve 
is shown. This half is obtained from introducing the branch cuts which connect 
the branch points —00 and 0, and s\ and (bold lines along the real axis). The 
signs below and above these branch cuts indicate the value of w 'right above' 
and 'right below' the branch cut. At '+', w has the value +i|w|; at ' — ', w has 
the value — i|w|. The integration path 7^ is a closed loop which encircles the 
right branch cut. The integration path 7^ consists of two parts. It starts at 
z = 1 on the shown Riemann sheet of w and intersects the branch cut between 
and —00. This part is marked by a solid line. The integration path continues 
on the other Riemann sheet (which is a copy of the one shown and which joins 
this copy at the branch cuts) and ends at the point 1 on the other sheet. This 
part is shown as the dashed line, (c) The continuation of (b) for s\ leaving the 
real axis (see Sec. H]). 
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Figure 10: Effective potentials and energies (a) and phase portraits (b) for 
constants of motions (si,S2) (or equivalently {k,l)) in the regions BBi, BB2, 
BB3, WGi, WG2 and WG3 defined in Fig. [6l For the ^ degree of freedom, the 
hatched regions mark the forbidden regions [0,/iB] and [2K{q') — fiB,'2K{q')] 
which are not contained in the region ([3]). ((a^,c^) = (3/2, 1/2).) 



with the closed integration path 7^ in (|56p defined as in Fig.[3^b). This assigns 
a finite, positive real valued integral also to the unbounded degree of freedom 
A or equivalently ^ for motion T2. For motion T3 the order of S2 and along 
the real axis in Fig. [9ljb) is reversed. The integral defined according to ((56|) 
is then negative real. Though at first not important for the classical dynamics, 
this integral will play an important role in the semiclassical computations in 
Sections [5] and [HI 



4.2 The 3D system 
4.2.1 Phase space foliation 

Similarly to (P7|) the separated momenta conjugate to (C,'7jC) can be written 
for the 3D system as 

where (s g {C, ?7, £,}) (see 52): ^^'^ ' separation constants, and < S2 are 
defined as in . Since the latter are the squares of the zeroes of the numerator 
polynomial on the right of (|58p . they are the squares of the turning points in 
the respective degree of freedom s G {-q, (, ^}. 
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The corresponding equations for the coordinates {i^, ^, A) are 



pl = as^is\s)~2ks\§) + l), (59) 

where s G {1^,11, X}, and s(s) e {(^{h'),r]{p), ^(X)} are the functions defined in 
([25]) . Similarly to the 2D case, ([551) and ([5^ are the classical analogues of the 
separated wave equations ([^^ and ([55]) , respectively. The specular reflection at 
the hyperboloidal boundary ry = f or equivalently 11 = fiB and fi = 2K(q') — /iB 
becomes 

iC,v,£„PC'Pv^Pi) '-^ {C,v,£„Pc^-Pv^Pi) (60) 

or 

{i^,fi,X,p^,p^,Px) {iy,fi,X,p^,-p^,px), (61) 

respectively. Note that, apart from the specular reflection, the motion described 
in terms of the phase space coordinates {s,ps), s € {v, ^, X} is smooth on the 
double cover ((5^ . 

Expressing the separation constants k and I, or their energy scaled coun- 
terparts K := 2mEk and L := 2mEl, in terms of Cartesian coordinates and 
momenta gives 

K = li\L\' + {a' + c')pl + aV. + c^pl), (62) 
L = c^Ll + a^Ll + a^c^pl, (63) 

where , Ly and denote the components of the angular momentum about 
the origin L = r x p. The separation constants together with the total energy 
^ = {Px + Py +Pz)/(2"i) give three constants of motion. Hence, the classical 
system is integrable, and a modification of the Liouville- Arnold theorem implies 
that the six dimensional phase space is foliated by the common level sets of E, 
K and L, or equivalently i?, s\ and which are invariant cylinders. 

Like in the 2D case the energy plays no major role since it only determines the 
speed of the motion along the straight lines (in configuration space) . So in order 
to describe the foliation of the energy surfaces by the invariant cylinders it is 
sufficient to consider the energy surface of a fixed positive energy E > Q. Other 
foliations of other energy surfaces are then obtained from a suitable scaling. 
On such an energy surface, there are then smooth two parameter families of 
invariant cylinders parametrized by s\ and s^. The parameterization intervals 
of these cylinders can be obtained from requiring the momenta ()58p to be real 
and analyzing the disposition of zeroes s\ and s| relative to the poles a? and 

in (j58p . To obtain real momenta, s\ and hence s\ can only take nonnegative 
values. Similar to the 2D case we will therefore occasionally use si = \fs\ and 
S2 = \/sf. It then turns out that there are six different smooth families of 
invariant cylinders which we denote by BBi, BB2, BB3, WGi, WG2 and WG3 
as shown in the bifurcation diagram in Figures [61(b) and (c). 

In order to describe the motions on the different families of cylinders it is 
again useful to illustrate the corresponding effective energies and potential, and 
phase portraits and also the intersections of these cylinders with the various 
Cartesian coordinate planes. This is shown in Figures fTUl and [TT| respectively. 
To simplify the discussion we will consider, like in the 2D case, the specular 
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reflections at the billiard boundary to be smooth. In the 3D case this implies 
that we identify and — when /x = /iB or /i = 2K{q') ~ /ib (see Fig. [TUlb)). 

For a fixed energy E > a pair {si, S2) (or the corresponding pair (A:, I)) in 
BBi or BB2 has as its level set a toroidal cylinder R x which we illustrate 
in terms of its projection to configuration space in Fig. 1111 It is unbound in 
the direction of A and the motion is oscillatory in the transverse degrees of 
freedom 1/ and n. In BB2 the motion oscillates with reflections at the boundary 
hyperboloid. The intersection of the cylinders of type BB2 with the y—z plane is 
bounded by the two branches of the hyperbolas rj = si, similar to the "bouncing 
ball modes" which one flnds in the billiard in a planar ellipse [23]. 

In contrast to that the motion in BBi, though oscillatory in 77 and does not 
touch the boundary hyperboloid, i.e. the corresponding toroidal cylinders are 
foliated by straight lines of free motions without reflections. A pair (5^,32) in 
BB3 represents motion which does not cross the y—z plane. The corresponding 
level sets consist of two toroidal cylinders R x which are bounded away from 
the y — z plane by the ellipsoid £, — 82- 

Pairs (s^jS^) in WGi or WG2 involve motions which are rotational in ( 
(or, equivalently, in v). They represent two toroidal cylinders R x which 
differ by the sense of rotation (see the corresponding panels in Fig [TT|) . In the 
elliptical cross-section in the x — y plane the motion WG2 is bounded by the 
ellipse 77 = si, similar to the "whispering gallery modes" which one finds in 
planar elliptic billiards. As in the case of BBi, motions in WGi do not touch 
the hyperboloidal boundary. The corresponding toroidal cylinders are again 
foliated by lines of free motion without reflections. For (s^,S2) in WG3 the 
rotational motions are again bound away from the y — z plane by the ellipsoid 
^ = S2- The corresponding level set consists of four toroidal cylinders which 
have a; > or X < combined with different senses of rotation. 

The smooth families of cylinders bifurcate along the boundaries of the s\ — 
s\ parameterization intervals in Fig. [H] Along S2 = 1 we have the (minor) 
bifurcation from cylinders consisting reactive trajectories which have reflections 
at the boundary hyperboloid to cylinders with trajectories having no reflections. 

Along s\ = (? the motions bifurcate from bouncing ball to whispering gallery 
type. Note that the distinction between BBi and WGi is only 'artificial' . They 
both consists of free motions (without reflections). For such motions, there are 
more constants of motion than degrees of freedom (the free motion can be sepa- 
rated in several coordinate systems). In such so called superintegrable systems 
(a multidimensional harmonic oscillator is a simple example) the foliation by in- 
variant cylinders (or equivalently invariant tori in the case of a compact system 
like a harmonic oscillator), is therefore not uniquely defined ^9J. 

Most importantly for the reaction dynamics is the bifurcation from nonreac- 
tive motions to reactive motions along s\ = a^. In fact, the joint level set of the 
two constants of motion E and S2 for a fixed E > and S2 = consists of the 
energy surface of the unstable invariant two-degree-of-freedom subsystem which 
consists of the billiard in the bottleneck ellipse which has x — and Px — and 
its stable and unstable manifolds. We illustrate the foliation of this level set 
in Fig. [T21 The two-degree-of-freedom billiard in the bottleneck ellipse can be 
viewed to form the transition state for the 3D system. The transition state at 
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energy E is then given by 

TS = {{x,y,z,px,py,p;,) e M.^ : X ^ 0, ^ 0, 

y^ + ^ <l,pl+pl^27nE}. '^^^^ 

The biUiard in an elhpse is foUated by two different smooth famihes of twodi- 
mensional tori which in this case are parameterized by s^. For < sf < c^, 
these tori are of bouncing ball type, and for < < 1, the tori are of whisper- 
ing gallery type. The bifurcation between theses two families at = involves 
the unstable periodic orbit along the major axis of the bottleneck ellipse. At 
— the bouncing ball tori degenerate to the stable periodic orbit along the 
minor axis of the ellipse. At = 1 the whispering gallery motions degenerate 
to the two periodic orbits sliding along the perimeter of the ellipse in opposite 
directions (see [23] for a detailed discussion). Regarding the specular reflections 
to be smooth, the energy surface of this invariant subsystem forms a three- 
dimensional sphere, SS^. In the full original 3D system this sphere is unstable 
with respect to the transverse directions parametrized by A and px and therefore 
has stable and unstable manifolds and W" which are also contained in the 
level set (see the lines in the X — p\ plane in the left panel of Fig. fT^ . The 
topology of and can be inferred from taking the Cartesian product of 
the 3-sphere of the invariant subsystem with the lines in the X — p\ plane in 
the left panel of Fig. [l^l i.e. the stable and unstable manifolds have topology 
K, X SS"^. Like in the case of the 2D system discussed in Sec. 14.11 these stable 
and unstable manifolds are again of codimension 1 in the energy surface. This 
way they again have sufhcient dimensionality to act as separatrices, and in fact 
they again separate the reactive trajectories from the nonreactive trajectories. 
Similar to the case of the 2D system described in Sec. 14. li the manifolds W and 
VK" again each have two branches. We again denote the branch of which 
has Pa > (resp. p\ < 0) the forward (resp. backward) branch, (resp. W^), 
of the stable manifold. Similarly we again denote the branch of W'" which has 
Px > (resp. Px < 0) the forward (resp. backward) branch, (resp. W^) of 
the unstable manifold (see Fig. fT^ . Also, we again call the union of the forward 
banches, 

Wf := U Wf , (65) 
the forward reactive cylinder, and the union of the backward branches, 

Wb := U , (66) 

the backward reactive cylinder. These forward and backward reactive cylin- 
ders then again enclose the forward and reactive trajectories, respectively, and 
separate them from the nonreactive trajectories in the energy surface under 
consideration. 

Moreover, we can define a dividing surface DS by setting x = on the energy 
surface which gives 

^2 

DS = {{x,y,z,p^,py,p^) G : a; = 0, + -rj < 1, ^^^^ 

pI + pI + pI = 2mE} 

With our convention to consider the specular reflections to be smooth the di- 
viding surface DS has the topology of a four-dimensional sphere, S"^. Similar to 
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the situation in the 2D system the three-dimensional sphere associated with the 
transition state TS in can again be viewed to form the equator of the DS 
4-sphere. In fact the transition state TS divides the dividing surface into two 
hemispheres, the forward dividing surface 

DSf = {{x,v,z,v^,Vv,Vz) e R''^ X ]R3 : a; ^ 0, y2 + — < 1, ^^^^ 
pl+p^+p2^2mi;,p, >0}. 
and the backward dividing surface 



DSb = {{x,y,z,p^,py,p,) e X ]R3 . 2; ^ 0, + ^ < 1, 



£2 - ' (69) 



pi+pt.+pi^2mE,p, <0}. 



Each forward reactive trajectory has a single intersection with the forward hemi- 
sphere; each backward reactive trajectory has a single intersection with the back- 
ward hemisphere. Nonreactive trajectories do not intersect the dividing surface 
DS at all. Like in the 2D case the dividing surface is everywhere transverse to 
the Hamiltonian flow apart from its equator which is an invariant manifold. 

4.2.2 Action Integrals 

As we have seen in the previous subsection the phase space of the 3D system 
is foliated by six different families of invariant toroidal cylincers, R x T^. For 
the toroidal base, which is associated with the degrees of freedom ( and 77 (or 
equivalently v and fi), we can again define action-angle variables. The actions 
in this case are given by 





- 2fcs2 


+ 1 


(s2- 


a2)(s2 


-C2)' 


'(s2- 


S2)(s2 


-si) 


(s2- 


a2)(s2 


-C2)' 



(70) 



with s e {CiV} s-iid the ps being taken from (|58p . The integers rris and the 
integration boundaries s_ and s+ can be found in Tab. [T] 

For the actions of the symmetry reduced system, which we again denote by 
/s, we always have rus — 2, s € {C, 77}, i.e. 



■1 



is = V2mE- / psds . (71) 



Substituting z = ,s2 in Eq. ((701) shows that the action integrals /,;; and Z,, are 
both of the form 

L^V2^^ r (^-^^)(^ -^i) d., (72) 
Att J^_ w{z) 

where 

5 

w^{z) = P5iz):^l[{z-z,), (73) 

i=l 
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Figure 11: Configuration space projections of the invariant cylinders corre- 
sponding to the motions BBi, BB2, BB3, WGi, WG2 and WG3 as their inter- 
sections (shaded regions) with the Cartesian coordinate plane. The bold lines 
mark the intersections of the boundary hyperboloid ([4]). The ranges for x, y 
andz are [-3/2,3/2]. {{a^,c^) = (3/2,1/2).) 
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Figure 12: Foliation of the level set E = const. > 0, S2 — which consists 
of the invariant billiard in the bottleneck ellipse that forms the transition state 
of the 3D system, and its stable and unstable manifolds with branches W^^/f, 
and WJ,f^, respectively. In the X — px phase plane the invariant billiard appears 
as the dot at the origin. The corresponding phase curves in the ^ — plane 
and V — pi, plane are parametrized by s\ E [0, 1]. For < < c^, the phase 
curves are of bouncing ball type (green dashed curves); for < < 1, the 
phase curves are of whispering gallery type (blue dashed curves) . The separatrix 
between bouncing ball and whispering gallery motions has = (red dashed 
curves in the right panels). Note that the two pieces of the phase space curve 
corresponding to the bouncing ball motion are mapped onto each other by the 
involution 
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Table 1: Integration boundaries, s_ and s^, and multipliers in ()70|1 for the 
six types of 3D motion BBi, BB2, BB3, WGi, WG2 and WG3. 
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and z_ and z+ are consecutive elements of the set {zi = 0, Z2 — s^, ~ b^, Z4 = 
S2, Z5 = a^jZh = !}• Again Zb = 1 corresponds to the boundary hyperboloid. 
The differential dz/w{z) has the six critical points {zi, Z2, 2:3, Z4, Z5, cxd} which 
means that the integrals (Tf^ are hyperelliptic. There do not exist tabulated 
standard forms for these integrals like for the elliptic integrals ([52]) in the case 
of the 2D system. However we can again view them as Abelian integrals, which 
in this case are defined on the hyperelliptic curve 

r^={{z,w)eC^ : = P5{z)}, (74) 
which is an algebraic curve of genus 2. The integrals (|72p then become 

/. = V2^'^ / (z - sl){z - sD- , (75) 

where s G {(, rj} ^-nd the integration paths 7^ and 7^ for motions of type WG2 
and BB2 are shown in Fig. [131 For motions of type WG3 and BB3, the order 
of S2 and along the real axis in Fig.[T3]is reversed which does not affect the 
definitions of 7^; and 7,,. The integration path 7^ is a closed path on F^, and 
hence the integral Iq is a complete hyperelliptic integral. Due to the billiard 
boundary the integration path 7,, is not closed; the integral is an incomplete 
hyperelliptic integral. 

Similarly to (15^ , we can also define a closed hyperelliptic integral associated 
with ^, 

= j [z- sl)[z - sD- , (76) 

27r J^^ w 

where 7^ for motions of type WG2 and BB2 is defined in Fig. [T3] and leads to 
a real positive for these motions. For motions of type WG3 and BB3 the 
order of S2 and along the real axis in Fig. [T3]is reversed, and becomes real 
negative. 

Moreover we define the integral 

Jc„ = W2mE- / (z-zi)(z-Z2)— , (77) 

Ac. ^ 

where 7,^^ for WG2 and BB2 is also defined in Fig. [T31 The change of the order 
of S2 and again does not affect the definition of 7;;^. This way we get a real 
positive I(;n for WG2/3 and a real negative for BB2/3. 

The integrals and ^ will play an important role in Sections [5] and [6l 

5 Computation of the classical and quantum trans- 
mission from transition state theory 

In this section we compute the transmission probabilities for the classical and 
quantum transport from the regions x <ti —1 (the 'reactants' region) to the 
region x 1 (the 'products' region) in the geometries (P) (2D) and ^ (3D). 
In the quantum case we are interested in the cumulative reaction probability 
which is defined as 

N{E) = Trf{E)f^{E) , (78) 
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Figure 13: (a) The graph of the polynomial Psiz) defined in ([74|) for real sepa- 
ration constants sf and s"^ satisfying < < 1 < < (motion type WG2). 
(b), (c), (d), (c) Complex planes with definitions of the integration paths 7,^, 7^, 
7j and 7^,,, along which , and I(^ri in (|75p . ((7^ . and ((77)) are computed; 

(b) , (c) are for real sf and (scattering states of type WG2 and BB2 , respec- 
tively), and (d), (e) are the corresponding continuations of (b) and (c) when si 
and s| leave the real axis (resonance states). Similarly to Fig. [SI the complex 
planes can be viewed as one half of the hyperelliptic curve r^, (the Riemann 
sheet of one 'sign' of the square root w). To make the square root w well defined 
in the complex plane three branch cuts connecting consecutive branch points of 
w are introduced (bold lines). The left cut connects the branch points and 
00; the middle cut connects and sf; the right one connects s| and a^. The 
integration paths 7^,, and 7^ are closed loops which encircle the middle and 
right branch cuts, respectively. The integration path 7^ is also a closed loop 
of Fu,. In the picture shown it consists of two parts. It starts at either z — 
(WG) or z = (BB) on the other Riemann sheet of w (which is a copy of the 
one shown and which joins the copy shown at the branch cuts) and intersects 
the branch cut between and 00. This part is marked by a dashed line. The 
integration path then continues on the shown Riemann sheet where it ends at 
either z = (WG) or z = sf (BB). The integration path 7^ also has two parts. 
It starts at z = 1 on the other Riemann sheet of w and intersects the branch 
cut between and . This part is marked by a dashed line. The integration 
path continues on the shown Riemann sheet where it ends at the point 1. Like 
in Fig. [5]the signs -f or — indicate whether w = -|-i|w| or w = —11^1 'just above' 
or 'just below' the respective branch cut. 
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where T{E) is the transmission block of the scattering matrix at energy E (for 
references in the chemistry hterature see, e.g., [30]; in the context of balhstic 
electron transport problems (j78p is known as the Landauer-Biittiker formula 

[snEaiss]). 

According to its definition, N{E) can be computed from the scattering ma- 
trix. However, this is a very inefficient (and for the systems with many de- 
grees of freedom even infeasible) procedure since one has to determine all the 
state-to-state reactivities while N{E) is merely a sum over these reactivities and 
hence no longer contains information about the individual state-to-state reactiv- 
ities. Much effort has been and still is put into finding a computationally cheap 
method to compute N{E). In the chemistry hterature (see, e.g., [SU [35l [30] ) a 
method has been developed to compute N{E) on the basis of transition state 
theory where the classical transmission probability is computed from the flux 
through a dividing surface which for a given energy E separates the energy sur- 
face into a reactants and a products region. For such a dividing surface DS the 
flux from reactants to products can be computed as 

f{E)^ f [ S{E-H{q,p))F{q,p)Pr{q,p)dfqdh. (79) 

Describing the dividing surface by a zero level set of a function s on phase 
space, or more precisely a function which is negative on the reactants side of 
the dividing surface and positive on the products side of the dividing surface, 
F in (|79p is defined as the following composition of functions 

F{q,p) = ^ e o s o ^'Hi<l,P)\t=o = S«1^P)){^' H}iq,p) . (80) 

Here O denotes the Heaviside step function, is the Hamiltonian flow gen- 
erated by H acting for the time t, and {•, •} denotes the Poisson bracket. The 
function Pr in ([7^ is defined as 

P,(g,p) = lim eoso$*,(g,p), (81) 

t — >oo 

which acts as a characteristic function on the dividing surface. In fact, by 
construction Pr{q,p) = 1 if the trajectory through the point {q,p) proceeds for 
t oo to products which is the region where the function s is positive and 
Pr{q,p) = otherwise (see [8 for a detailed discussion). 
The quantum analogue of (fT^]) is given by 

N{E) ^ 27rhTT 5{E - H)FPr , (82) 

where 

F = -i[e^,i?] (83) 

and 

A- = lim e^^'eTse"*^*. (84) 

t — >oo 

Here Q o s denotes a quantization of the classical function 8 o s [51 . 

Like its classical analogue ([79]) the evaluation of ([82]) involves a computation- 
ally expensive time integration which is manifested in (|8ip and (|84p , respectively. 
The computational advantage of the transition state theoretical formulation of 
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N{E) over the original definition in ([75]) is therefore not obvious. In practice 
one can carry out the time integration only to a finite time. This time has to 
be large enough so that it can be decided that after this time the resulting tra- 
jectory (classical) or wavefunction (quantum) will stay in the products region. 
In order to minimize this integration time one has to choose a "good" dividing 
surface. In fact, for a dividing surface that, classically, is crossed exactly once 
by all reactive trajectories and not crossed at all by nonreactive trajectories (see 
our discussion in Sec. \^ no time integration is required at all. The character- 
estic function in ()8ip can then be replaced by a function which at a point 
{q,p) on the dividing surface is one if the Hamiltonian vector field at this point 
pierces the dividing surface in the forward direction and zero if the Hamiltonian 
vector field at that point pierces the dividing surface in the backward direction. 
In other words this means that we can omit the function P,. in (|82p and re- 
strict the integral ([5^ to the foward hemisphere of the dividing surfaces that 
we constructed in Sec. HI The choice of a good dividing surface is thus crucial 
to benefit from the transition state theoretical approach to compute classical 
and quantum transmission probabilities. 

In Sec. m we used the separability of the transmission problem discussed in 
this paper to construct the dividing surface which has the desired properties. As 
a consequence of this separability we similarly get in the quantum mechanical 
case that the transmission subblock of the scattering matrix in ([78| is diagonal. 
Using 

tnrn{E) — tn{E)6nm , (85) 

where n and m label the scattering states at the energy E, and Snm is the 
Kronecker symbol we have 

7V(i?) = ^r„(£;), (86) 

n 

with the state-to-state transmission probabilities defined as Tn{E) := \tnn{E)\'^ . 
In the following we present the computation of these transmission probabilites 
and the comparison of the resulting cumulative reaction probabilities with the 
classical fiux. 

5.1 The 2D system 

5.1.1 The quantum transmission 

To compute the transmission probabilities T„ in the 2D case we look for solutions 
of the form 

^C;n (0 V'? ;n (0 + r„^C;« (0%,n (0 (87) 

at the bottom {x <^ —1) and 

tnt^cACW^-niO (88) 

at the top {x ^ 1). Such solutions can be computed from first solving the 
transversal component C of the wave equations (jlip for the corresponding bound- 
ary conditions with £' as a parameter. As discussed in Sec. [3] the boundary 
conditions for -0^ at C = is determined by the paritiy tTj,. For TTy ~ +, we 
have ip'^{0) = (and choose ■(/'c(O) = 1); ^^'^ for tTj, = — , we have "0^(0) — 
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(and choose V'^(O) = !)• The Dirichlct boundary condition at C = 1 requires 
?/'^(l) — 0. This way we obtain modes which we label by the Dirac "kets" 
\ni^]TTy) where is a non- negative quantum number which gives the number 
of nodes of V'c in the open interval < C < 1- The modes of energy E deter- 
mine the separation constants Sj^^^.^ )(^)- This separation constant can then 
be used in the equation for the ^ component of the separated wave equations 
(fTTj) to find solutions of the form ([57]) and ([55]) . This however is not completely 
straightforward and does not give much insight into the structure of the solu- 
tions. We therefore resort to a semiclassical computation which will also lead 
to the semiclassical computation of resonances as we will discuss in Sec. |6l The 
semiclassical approximation is obtained from using S2(n^-7r )(-^) compute the 
transmission probability r(„^.^^)(£') as 

(a)T(„^.,^) (E) = — -i — — , (89) 

l-Hexp(6'(„^.^j(£;)/?i) 

where 6^nc^;-Ky) is a tunnel integral. This tunnel integral describes the quantum 
mechanical tunneling through the dynamical barrier which in terms of the A 
coordinate occurs as the barrier in the associated effective potential Vx^cff for 
motions of type T3 (see Fig. [7]). This tunnel integral is given by 




9(n,-n,){E)^-2iJ^ PAdA = -4iV2^/ J^^^^d^, (90) 

where S2 = ^2(n -tt )(-^) ^'^d A+ and A_ = — A+ are the corresponding turning 
points using the phase space coordinates {\,px) (sec, e.g., [36 for a derivation 
of the expression (|89p). For S2 > or equivalently E\ cff < V\^cs{0) which 
corresponds to classical reflection of type T3, px is imaginary along the inte- 
gration interval which is bounded by the real classical turning points A_ and 
A+ = — A_. This integral can be identified with two-times the integral of 
from a to the corresponding turning point S2 which gives the second equality in 
(^0)1 . For s| < or equivalently Ex,cS > Vx,cs{0) which corresponds to clas- 
sical transmission of types T2 and T3, the classical turning points X± become 
imaginary (with A_ being complex conjugate to A_|_) whereas px is real on the 
imaginary axis between X± . The branches of the square root in (|90p are chosen 
such that the tunnel integral is positive if < s\ and negative if s\ < a? . This 
choice of the branches can be described more precisely from relating ^jn^.^^-j (_E) 
to the integral that we defined in ([5^ . In fact we have, 

h^-^e(n,;.,){E). (91) 

The boundary value problem for -01^ can be solved numerically using a shoot- 
ing method which relates the solution of the boundary value problem to a New- 
ton procedure (see [37], and [131 [22] for similar apphcations). Since t^c is an 
oscillatory function which leads to multiple zeroes in the resulting Newton pro- 
cedure the shooting method requires good starting values s|. These are obtained 
from a semiclassical approximation also of the boundary value problem for . 
To this end we note that the phase portraits of the motions of type T2 and T3 
between which the classical motion switches from transmission to reflection are 
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identical in the v — pi, plane (see Fig. [7])(b). For these types of motions we can 
thus use the EBK quantization condition for the action Iq, 

= h{n^ + 1) , e Mo , (92) 

which is the same as the EBK quantization of a one-dimensional square well 
problem. We can also rewrite this quantization condition in terms of the action 
= of the symmetry reduced system which gives 

= h{n^ + i(3 - -Ky)) , e Mo . (93) 

This decomposes the semiclassical modes in terms of the parity -Ky. The quan- 
tum numbers uq and are related by 

n^ = 2n„ + i(l-7ry). (94) 

We note that for the type Ti the motions involve a smooth rather than a 
hard wall reflection in the v degree of freedom. As a result the EBK quantiza- 
tion for Ti would be different from the EBK quantization for T2 and T3, and 
hence, in order to describe the transission from T2 to Ti a unifrom semiclassical 
quantization scheme would be desirable. However this transition plays no role 
for the transition from transmission to reflection (see below) and we therefore 
do not consider this aspect in more detail. The quantization condition ()92p can 
be solved by a standard Newton procedure. The solutions for E and s| for a 
given quantum number tiq and parity tTj, are then used as the starting value for 
the shooting method described above. 

The cumulative reaction probability N{E) is then the sum over all the 
T(^n^;Try){E) in ()89p for all quantum numbers and parities TTy. For the nu- 
merical computation of N{E) we need only consider the finite number of modes 
which, at a value E > 0, have a nonnegleglibile transmission probability. A 
graph of N{E) is shown in Fig. [TH We note that on the scale of the picture 
one can notice no difference between the exact and the semiclassically computed 
N{E). Depending on the shape parameter for the boundary hyperbola the 
cumulative reaction probability shows more or less pronounced steps with unit 
step size. A detailed analysis of the graphs of N{E) can be obtained from re- 
lating the modes |ni^;7rj,) to the classical motions. For a given energy E, this 
relationship is established via the separation constant s'^^^^^.^ )(-^) which deter- 
mines the classical invariant cylinder the mode is associated with. As can be 
seen from the projections of the cylinders in Fig. [S] these projections become 
increasingly confined in the order T3 ^ T2 — > Ti. Since high confinement in 
configuration space implies high kinetic energy via the Heisenberg uncertainty 
principle, the modes, which classically correspond to the type of motion Ti have 
highest energy. In fact, for low energies all modes have S2(n^-7r )(-^) clas- 
sically reflecting type of motion T3. Upon increasing the energy the S2(-„^.^ ^ (E) 
wander towards the transmitting mode T2, and for even higher energy to Ti, 
see Fig. inja). Concerning the classical mechanics, the border between reflection 
and transmission is given by s| — . This border is crossed for the modes 
\nc^\ TTy) for different energies. Upon crossing the border the tunnel integral 1^0]) 
changes sign and the transmission probability ([5^ changes from to 1. The 
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Figure 14: Top Panel: Cumulative reaction probability N{E) as a function of 
the wavenumber k = \/2mE /h for the shape parameter a? — 5, and, for com- 
parison, = oo, which corresponds to the transmission through a rectangular 
strip. The ticks on the wavenumber axis mark the energies at which, for — 5, 
the modes |n^;7ry) "open" as transmission channels (see text) (the key to the 
tick labels is given in the table). The smooth dot-dashed blue curve and the 
solid blue curve show the Weyl approximations of N{E) defined in and 
([loop , respectively. Bottom Panel: Resonances in the complex wavenumber (k) 
plane, for — 5. Semiclassical resonances are marked by pluses (-1-) and exact 
resonances by diamonds (0). Note that we use the symbol k for wavenumber to 
distinguish it from the separation constant k of the 3D system, [h = 1, m — 1.) 



36 



energy for which the tunnel integral of a given mode |n^; tt^) is zero, and hence 
gives T(„^.^^) {E) = 1/2, can be defined as the energy at which the mode "opens" 
as a transmission channel. Marking these energies on the energy axis in Fig. 
we see in which order the transmission channels open and this way contribute 
a step of N{E). Semiclassically these "opening" energies are identical to the 
eigenenergies of a square well. 

5.1.2 The classical transmission 

The classical transmission probability can be computed from the directional 
flux through the dividing surface DS of energy E defined in (|46p. or following 
our discussion at the beginning of this section by an integral over the forward 
hemsiphere DSf of this dividing surface. In a more modern notation which also 
reveals the symplectic nature of the flux (see [3H1[3H] and also f2Dl) the flux is 
given by 

f{E)^ f Lu, (95) 

JDSf 

where lu is the symplectic 2-form 

uj = dx A dpx + dy A dpy . (96) 

Since lu = d(j) where (f) is the Liouville 1-form 

(f> = Px dx + py dy (97) 

we can utilize Stokes' theorem to compute f{E) from integrating cj) over the 
boundary of the forward hemisphere DSf. Using the fact that the boundary of 
DSf is given by the transition state TS consisting of the periodic orbit along the 
y axis at energy E (see Sec. I4.ip we find that the flux is given by the Liouville 
action of the periodic orbit, 

f{E) = [ (f) = 4V2to£; . (98) 
Jts 

In order to make the comparison to the cumulative reaction probability N{E) 
we consider the dimensionless quantity 

NwcyiiE) = ^f{E) , (99) 

which is shown together with N{E) in Fig. [T31 We see that A'wcyK-E) gives 
an approximate smooth local average of N{E) which however overestimates the 
local average of N{E) as the graph of A'wcyi(-E') intersects the graph of N{E) 
at the top of its steps. In fact disregarding the tunneling, N{E) simply gives 
the integrated density of states of the transition state or activated complex (the 
one-dimensional square well along the y axis) to energy E. The term Ny/(.y\{E) 
is the Weyl approximation of this quantity. As also shown in Fig. [14] one can 
obtain a better local average by modifying N^eyi{E) to 

NwcyiiE) = ^fiE) - \ , (100) 

which can be formally derived from counting the mean number of states to 
energy i? in a one-dimensional square well potential. 



37 



5.2 The 3D system 

5.2.1 The quantum transmission 

To compute the transmission probabilities T„ in the 3D case we look for solutions 
of the form 

^C;n(C)V'r,;n(?7)V'«;n(C) +^«V'C;«(C)V'r,;n(??)V4„(e) (101) 

at the bottom (x ^ — 1) and 

tn^(AC)MvW^;niO (102) 

at the top (x ^ 1). In this case we first solve the components of the sepa- 
rated wave equations ()24|) and the corresponding boundary conditions which 
belong to the transversal coordinates C and rj with the energy i? as a param- 
eter. The boundary conditions for ipr/ sue given by the parity tt^ which yields 
the index oi ijjjj at ij — c and the Dirichlet boundary condition — 0. 

The boundary conditions for ipi^ are determined by the parities tt^ and tTj^: tt^ 
determines the index of ipi^ at C = c and tTj, determines whether ^/'^(O) — 0, 
V'c(O) = I {TTy = +) or Vc(0) = 0, V'c(O) = 1 (tTj, = -)• This way we 
obtain modes that are parametrized by E and which we label by the Dirac 
kets n^; TTj,, TTz), where n,^ and are non-negative quantum numbers which 
give the number of nodes of ip^ and ip,^ in the open intervals c < 77 < 1 and 
< C < c, respectively. The modes for energy E determine the separation 
constants (fc(„^,„,j;^^,^^)(£;), /(„^,„^;^^_^^)(i;)). These can then be u sed in the ^ 
component of the equations to find solutions of the form (jlOip and (|102p . 
Like in the 2D case we resort to a semiclassical computation of the transmission 
probabilities instead. Analogously to we obtain 

T(„,„,;.„,..)(i?) = i + exp(0(„,,,'^.^,.,)(£;)/?.) ' 
where 0(^nf,nn:iry,Tr^)iE) IS the tunnel integral 

in,,n,;.,^.^){E) = -2i f %xdX (104) 

S2 



which describes the tunneling through the potential barrier of the effective po- 
tential V\^cS for types of motion WG2/3 and BB2/3 (see the corresponding phase 
portraits in Fig. [TU|) . The branches of the square root in (|104p are again chosen 
in such a way that the tunnel integral is positive if < s\ (corresponding to 
classical reflection of type WG3 or BB3) and negative if s\ < (corresponding 
to classical transmission of type WG2 or BB2). We can again make this more 
precise by relating 9(n^,nn\TTy.-K,){E) to the integral that we defined in (ffH]) . 
This gives 

k = -^(^in„n,:.y..AE). (105) 

Like in the case of the 2D system we solve the boundary value problems for 
and V't? by a shooting method. To this end we again need good starting 
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values for the separation constants l(n^,nn;Try.-K^)iEI) and fc(„j,n^;7r^.7rj which 
we obtain from a semiclassical computation. In contrast to the 2D system we 
here face the problem that the types of motion WG2/3 and BB2/3 differ with 
respect to their degrees of freedom ( and 77, or equivalently ly and /i (see the 
corresponding phase portraits in Fig. [T0|). Accordingly, the EBK quantizations 
of the actions are different. The conditions are 

Iu = {n,. + ^)h, I^ = {n^ + ^)h (106) 

for the bouncing ball motions of type BB2/3, and 

3 
I^ = {n^ + -)h, I^^{n^ + -)h (107) 

for the whispering gallery motions of type WG2/3. Writing these EBK quanti- 
zation conditions in terms of the actions of the symmetry reduced system one 
finds 

/C = (nc + i(2-7r,))ft, i^ = {n, + ^{3-7r,))h (108) 
for the bouncing ball motions of type BB2/3, and 

i^^{nc + ^i2-7ry-TT,))h, i^ = {n^ + ^i3))h (109) 

for the whispering gallery motions of type WG2/3. The quantum numbers 
(7i,y,7i^) of the full system and the quantum numbers (71^,71,,) of the symmetry 
reduced system are related by 

^ 2nc; + ^{2 ~ TTy) , = + ^(1 - tt^) , (110) 

for the bouncing ball motions and 

TT-t, = 27T.<; + i(2 - TTy - TT^) , 71^^=71,,, (111) 

for the whispering gallery motions. We can overcome this problem of differing 
quantization conditions by introducing a uniform quantization of the actions 
and Irj which interpolates the EBK quantizations in the regions BB2/3 and 
WG2/3 in a smooth way. Using ideas similar to [52] one finds that (jlOSp and 
(|109p can be written in the uniform way 

ic^inc + ^)h, i,^{n, + ^)h, (112) 

where again n^, G INq, and and are "effective Maslov indices" given by 

oir = TTz— arctane^'^^'' + 2 — TTy — TTj , (113) 

TT 

2 

a„ = TT^- arctane^"/'* + 3 - TT^ . (114) 

TT 

Here 9^ and O^i are again tunnel integrals which in this case describe the tun- 
neling through the barriers of the effective potentials V^^cS and V^^cS in Fig- UHl 
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that one needs to overcome to change between whispering gallery and bounc- 
ing ball type motions (see [53] for more details). The tunnel integrals 6*^ and 
Or, are again best described as integrals on the hyperelliptic curve in ((7^ . 
Interestingly 9q = —Or, and 

k = -S^^~^kn, (115) 

where 7,^^ is the integral we defined in ((77| . For Oq = —Or, ^ fi, and using (|110p . 
one recovers the quantization conditions for the bouncing ball type motions 
in (|106p : for —0(^^ 0,, S> and using piip . one recovers the quantization 
conditions for the whispering gallery type motions in (jl07p . 

In contrast to the other motions, we note that the types of motion BBi and 
WGi involve a smooth rather than a hard wall reflection in the /i degree of free- 
dom which leads to yet another set of EBK quantization conditions. However, 
as we will see below, for the energies under considerations BBi and WGi play no 
role for the transmission problem (see the discussion for the analogous effect in 
the 2D in Sec. 15.1.1]) . The uniform quantization conditions pi2p can be solved 
by a standard Newton procedure. The resulting values for s\ and s| for given 
quantum numbers uq and n^, and parities 'Ky and tt^, are then used as the start- 
ing values for the shooting method to solve the C and rj components of the wave 
equations as described above, and hence to compute the transmission probabil- 
ity in (|103p . The cumulative reaction probability N{E) is the sum over all these 
transmission probabilities. As in the 2D case, to numerically compute N{E) we 
need only consider the finite number of modes which, at a value _B > 0, have a 
nonnegleglibile transmission probability. A graph of N{E) is shown in Fig. [TSl 
Depending on the shape parameters (a^,c^) for the boundary hyperboloid the 
cumulative reaction probability shows more or less pronounced steps which in 
contrast to the 2D case (see Fig. fH)) are of size 1 or 2. This can be understood 
in more detail if we relate the modes jn^;, n^; tTj,, tt^) to the classical motions. For 
a given energy E, this relationship is established via the separation constants 
{k{n^,n^;Try,jr,}{E),l(n(,n^;7ry,TT,){E)) which determine the corresponding toroidal 
cylinders. The wave functions of the modes jriij, n,,; tTj,, tt^) are mainly "concen- 
trated" on the projections of the corresponding toroidal cylinders to configura- 
tion space. As can be seen in Fig. [Til for the whispering gallery types of motion, 
these projections become increasingly confined in the order WG3 WG2 
WGi. For the bouncing ball types of motion the confinement increases in the or- 
der BB3 BB2 BBi. Since high confinement in configuration space implies 
high kinetic energy via the Heisenberg uncertainty principle, the modes, which 
classically correspond to the types of motion WGi or BBi, have highest energy. 
In fact, for low energies all modes have (fc(„^,„„;7r„,7r,)(£'), ^(«c,n^;7r^,7r,)(^)) in 
the classically reflective types of motion WG3 or BB3. Upon increasing the 
energy the (fc(„^,„^;^^_^^) (7;), ^(„^,„^;7r^,7r,)(£')) wander towards the transmitting 
modes WG2 or BB2, and for even higher energy to WGi or BBi, see Fig. [6] 
Concerning the classical mechanics, the border between reflection and transmis- 
sion is given by S2 — or I — — (a^ — fc)^- This border is crossed for the 
modes n,,; tTj,, tt^) for different energies. Upon crossing the border the tun- 
nel integral changes sign and the transmission probability changes from to 1. 
The energy for which the tunnel integral of a given mode rt^; tTj,, tt^) is zero, 
and hence gives Ti^ri^,n^:-n:y.T!^){E) = 1/2, can be defined can be defined as the 
energy at which the mode opens as a transmission channel (see the analogous 



40 




Figure 15: Top Panel: Cumulative reaction probability N{E) as a function 
of the wavenumber k = y2rriE/h for shape parameters (a^,c^) = (5,0.2), 
and, for comparison, (a^,c^) = (oo,0.2), which corresponds to the transmission 
through a cylinder with elliptical cross-section. The ticks on the wavenumber 
axis mark the energies at which, for (a^,c^) = (5,0.2), the transmission chan- 
nels n,;; TTj,, TT^) "open" (see text) (the key to the tick labels is given in the 
table). For pairs of near degenerate states the one corresponding to the higher 
wavenumber is marked above the k axis. The smooth dot-dashed blue curve and 
the solid blue curve are the Weyl approximations of N{E) defined in p25p and 
(|126p . respectively. Bottom Panel: Resonances in the complex wavenumber (k) 
plane, for — 5. Semiclassical resonances are marked by pluses (-I-) and exact 
resonances by diamonds (0). Note that we use the symbol k for wavenumber 
to distinguish it from separation constant k. (h — I, m = I.) 
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Figure 16: The {k, ^)-spectra of the modes leading to the jumps of the cumu- 
lative reaction probabihty in Fig. [15] for (a^,c^) = (5,0.2). For each shown 
mode n^; TTj,, TT^) (see the table in Fig. I15p the energy E is varied from 2 (for 
which (fc(„^^„^.^^_^^) (£;),?( \{E)) are in the reflective types of motion 

BB3 or WG3 beyond the right border of the shown region) to i5 = 100 for 
which {k(n^^n^.^y,-R,){E),l( ){E)) are in one of the classically transmit- 

ting types of motion BBi, BB2, WGi or WG2. The bold lines mark the classical 
bifurcation diagram. 
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Figure 17: Probability contours in the section a; = of the wave functions of 
the modes 1?^^ , h^;; tt^, tt^) at the moment when they "open" as transition chan- 
nels (see text). The wave functions are displayed in the order they contribute 
to the cumulative reaction probability in Fig. 1151 Light blue corresponds to 
low probability, red corresponds to high probability. The black lines mark the 
envelopes of the corresponding classical motion. Ellipses indicate whispering 
gallery modes; hyperbolas indicate bouncing ball modes. 



definition for the 2D case). These energies are marked on the energy axis in 
Fig. [13 Semiclassically these opening energies are identical to the eigenenergies 
of the ellipse billiard. 

Classically, the border S2 — corresponds to the unstable invariant motion 
in the y — z plane. This is the planar billiard in the bottleneck ellipse which 
is an invariant subsystem with one degree of freedom less than the full three- 
dimensional billiard. 

Due to the dynamical barrier the wave functions of the modes deep in the 
reflective types of motion WG3 and BB3 have negligible amplitudes in the y — z 
plane. As the energy increases the increase of the amplitudes is indicated by 
the switching of the corresponding transmission probability T(^nc^.nr,:vy,-n,) from 
to 1, i.e. the "opening" of a new transmission channel. The wave functions of 
the transmission channels which lead to the step in Fig. [Tulare shown in Fig.fTTI 
as their intersection with the y — z plane. 

The quantum mechanical manifestation of the two senses of rotation in 
the whispering gallery types of motion is the energetic quasi-degeneracy of 
the corresponding modes jn^;, n^; tTj,, tt^). The further the separation constants 
(/s(nf,n^;7r^,77,)(-E^), ^(n^ ;7r„ ,7r, ) in the whispcriug gallery types of motion lie 
away from the border s\ = to bouncing ball motions, the higher the effec- 
tive energy E^^cS lies above the effective potential V^.off. In this limit the role 
of the potential becomes negligible and the energy is essentially determined by 
the total number of nodes of ilii, along a complete v-\oo\) which is an ellipse in 
Fig. [T71 The relation (jllip between the quantum numbers of the full system and 
the quantum numbers of the symmetry reduced system leads to the energetic 
degeneracy of the two pairs of modes 

\nQ + l,nn,+,+) ^ jn^;, n^, -, -) , (116) 
\nQ,n^,+,-) ^ . (117) 

In Fig. [m this effect is seen for the pairs of modes |1, 0, +, and |0, 0, — , — ), 
|l,0,-,+) and|l,0,+,-), |2,0,-h,-H) and |1, 0, -, -), |2,0,+,-) and |2, 0, +), 
|2,0, — ,— ) and |3, 0,+,+) which, on the energy axis, become more and more 
indistinguishable as energy increases, and this way effectively lead to steps of 
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size 2 (see also Fig. [TBI and the plot of the wavefunctions in Fig. \T7\ . There is 
no analogous degeneracy for the bouncing ball type states, as can be deduced 
from the relation (jllOp . 



5.2.2 The classical transmission 

In order to compute the directional flux through the fourdimensional dividing 
surface ([57|) of the 3D system we consider the symplectic 2-form 

U! = dpx Adx + dpy Ady + dpz A dz , (118) 

from which we can define the 4-form 

ri' = icj^ ^dpx^dxAdpyAdy + dpx^dxAdpzAdz + dpyAdyAdpzAdz. (119) 

The directional flux through the dividing surface ([B7|) is then obtained from 
integrating W over the forward hemisphere DSf defined in i.e. 

f{E)= [ n'. (120) 

JDSf 

Noting that 

n' = d(j), (121) 

where is the 3-form 

^ 2^Pxdx + pydy + pzdz) Au; (122) 

we can again use Stokes' theorem to compute the flux from the integral over 
the boundary of DSf which is the transition state TS consisting of the invariant 
billiard in the bottleneck ellipse at energy E. Hence, 



fiE) ^ <P. (123) 



Using either (fHO)) or (fT^ we get 



f{E) = 2m7rVl - c^E , (124) 



which is the product of the area, A = ttvI— <?, of the bottleneck ellipse and 
the area of the circular disk of radius \/2mE in the two-dimensional momentum 
space {py,Pz) e 

In order to relate the flux to the cumulative reaction probability we introduce 
the dimensionless quantity 

1 A 9m F 

Comparing the graphs of iVwcyi(-E) and N{E) in Fig. [15] we see that Nwcyi{E) 
overestimates the local average of N{E). This is an indication that quantum 
effects are quite severe in this system. Using the fact that, neglecting quan- 
tum mechanical tunneling through the dynamic barrier, N{E) is essentially the 
number of states of the billiard in the bottleneck ellipse to energy E we can 
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introduce correction terms to -/Vw'cyi(-E') of which the first is proportional to \J~E 
and depends on the length, L, of the perimeter of the boundary ellipse and the 
second is a constant term resulting from integrating the Gauss curvature along 
the perimeter of the bottleneck ellipse [H] . This way we get 



, , A 2mE L y/2mE 1 

+ (126) 

where L = 4E{c) with E{c) denoting Legendre's complete elliptic integral of the 
second kind with modulus c. The graph of iVwcyi(-E') is also shown in Fig. [TS] 
and in fact gives a very good local average of N{E). 



6 Quantum resonances 

In Sec. m we have seen that the classical systems possess invariant subsystems 
of one degree of freedom less than the full system contained in the respective 
phase space bottlenecks. These systems which form the transition states were 
given by a one-dimensional billiard in a square well in the 2D case and the 
invariant elliptic billiard in the 3D case. The Heisenberg uncertainty relation 
rules out the existence of analogous invariant subsystems in the corresponding 
quantum mechanical problems. In fact, a wavepacket initialised on such an 
invariant subsystem will decay exponentially fast in time. This exponential 
decay is described by the resonances [8] . 

The resonances can be formally defined as the poles of the meromorphic 
continuation of the transmission probabilities to the lower half of the complex 
energy plane. A semiclassical approximation of the resonances can thus be 
obtained from the poles of the expressions of the transmission probabilities we 
have given in ([55)1 (for the 2D case) and p03p (for the 3D case). This leads to 
the complex EBK type quantization condition for the tunnel integrals 9 defined 
in (IHOl) or PM| given by 

9^mn{2nx + l),nxe'SMo, (127) 

or equivalently 

I^ = -ih{nx + ^),nxelNo, (128) 

(see Equations (pi]) and (|105p .) In the following we study this semiclassical 
approach to the computation of resonances and compare it with a numerical 
computation of the exact resonances based on the complex scaling method. 



6.1 The 2D system 

6.1.1 Semiclassical computation of resonances 

We at first compute the resonances semiclassically. To this end we have to 
simultaneously solve the (standard) EBK quantization condition ([55]) . i.e. 

/c = ^K + ^(3-%)), ncGlNo, (129) 
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in combination with the complex EBK quantization condition (jl28p which rewrit- 
ten in terms of the action of the symmetry reduced system defined in (j57p 
gives 

4 = -ih{n^ + i(2 - TT,) , e Mo . (130) 

As we will see below when discussing the numerically exact resonances this 
quantization condition decomposes the resonance states with respect to their 
parity tt^ ■ Note that similar to the semiclassical computation of the cumulative 
reaction probability in Sec. lS.l.Tl we assume that the type of motion Ti also plays 
no role for the computation of the resonances. Otherwise (|129|) would have to 
be replaced by a uniform quantization condition which interpolates between Ti 
and T2, i.e. across = 1 (see Sec l4.ip . This assumption is justified by the 
fact that the resonances are associated with the activated complex consisting 
of the classically invariant billiard in the one-dimensional square well potential 
which has =0^1 ^nd the resonances can be expected to have values near 
and hence stay away from s| = 1 . We will see that this assumption is indeed 
fulfilled. 

The solutions E and S2 of the quantization conditions (|129p and (|130p are 
complex valued. The integration paths 7,;; and 7^ defining I,^ and in (j55p 
and ()56p therefore have to be continued accordingly into the complex plane (see 
Fig. [9])(c). They can be found numerically using a standard Newton procedure. 
To this end one has to decompose (|129p and (|130p with respect to their real and 
imaginary parts which leads to a four-dimensional Newton procedure. We note 
that this procedure is less robust than in the real case in Sec. lS.l.Tl In particular, 
the procedure struggles when the energies E are close to the imaginary axis. For 
fixed quantum number and parity tTj, we find the resonances by starting near 
to, but not at, the corresponding "opening" tick on the real energy axis of 
Fig. 1141 and by smoothly moving the parameter into the complex plane. We 
then go through the grid (nj, tTx) & INq x {— 1, 1}. We give a list of the resulting 
complex energies E and separation constants S2 in Tab. [H 
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Table 2: The exact resonances {Eqm,S2 qm) and the semiclassical resonances 
(i^sc, s|_g(,) of the 2D system for Rei? < 50 and {u^^tTx} = {0,±} (the first two 
famihes of resonances). The relative error AE — {\Esc \ — |-Eqm|)/|-Eqni| is given 
in percent, {h = 1, m = 1.) 



6.1.2 Exact computation of resonances 

We compute the (numerically) exact resonances from the complex scaling method 
[m H3] . The main idea here is to turn the wavefunctions which are associated 
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Figure 18: Real parts of the separated resonance wavefunctions (solid lines) of 
the four states |n^, n^; ttj,, ttx) with uq = 2, £ {0,1}, % = — , and ttx = ±. 
For fixed {nf,7ry} (and hence fixed n^), and different and tTx, the wavefunc- 
tions tp^ are qualitatively the same (upper left panel). For increasing n\, the 
four right hand panels show the increase in the number of nodes of the separated 
wavefunction in the reaction coordinate a. The dashed lines arc the real parts 
of the corresponding effective potentials Vs^es{s), s € {i^, A}. We add the real 
parts of the effective energies Eg^eS (dotted lines) to the amplitude of the wave- 
functions to visualize their energies relative the height of the potential barrier. 
The bottom left panel is the real part of a sample effective potentials Va,eS over 
a larger a interval. 



47 



with the resonances and exponentiaUy divergent when the reaction coordinate 
goes to infinity into square integrable functions by a complex scahng of the 
reaction coordinate. In our case the reaction coordinate is given by ^, or equiv- 
alently A. In order to apply the complex scaling method we at first transform 
the A component of the separated Helmholtz equation (see (fT7|) ). 

'^^'^^ ^E(coshX^-sl)i>^. (131) 



2m dX 



The goal of the transformation is to get a wave equation which no longer in- 
volves an exponentially decreasing potential which would cause problems in the 
complex scaling method. This can be achieved by rewriting (|13ip in terms of 

cr(A) asinhA (132) 

and the scaled wavefunction -tjjcr defined by 

i/dcr/dA (cr-^+a^)^/^ 

This gives 

with the effective energy and potential given by 

Ea,M = E, (135) 

= a2 + a2 -2^A{<j^ + a^Y^ ^^^^^ 

respectively. As opposed to the effective potential in (|13ip . the potential Va^cS 
goes to zero as |cr| — > oo, and accordingly, for real E and s\, the wave function 
Tpa is a plane wave as ct — > oo. In fact, in terms of the original Cartesian 
coordinates [x, y) the solutions have to become plane waves for \x\ — > oo. From 
([5]) and noting that ^ = a cosh(A) we see that a = a sinh(A) is proportional 
to X and this is the motivation for the transformation (|132p . The scaling of 
the wavefunction (|133p is performed in order to again obtain a system of type 
'kinetic-plus-potential'. 

In order to compute the resonances we substitute a in ()134|) with ere'". 
Upon this scaling the outgoing plane waves become asymptotically decreasing 
as cr — > oo provided that a > — arg(i?). 

We implement the complex scaling method numerically using a shooting 
method. To this end we choose a suitably large but finite value Uoo at which we 
require the scaled wave function ijja to vanish. The other boundary condition 
on V. is given by V<t(0) = 0, ^/.^(O) = 1 if tt, = -, and V^,(0) = 1, ^',{0) = if 
TT^, = +. The boundary conditions for are the same as for the scattering states 
in Sec. 15.1.11 In order to implement the complex scaling method we decompose 
the two equations into their real and imaginary parts. This then leads to a (real) 
four-dimensional Newton procedure on the complex two dimensional E — s\ 
plane. As the starting values we use the semiclassical values for the resonances 
obtained as described above in Sec. 16.1.11 We note that the complex scaling 
method is quite sensitive with respect to the choice of aoc and the scaling angle 
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a. Like in the semiclassical computation this is particularly true for values of 
the energy near the imaginary axis. For the scaling used in these systems a 
"suitable" value for CToo ranges from « 1 for resonances with Re(-E) « 500 to 
« 5 for resonances with Re{E) « 3. An inappropriate value leads to apparent 
but false convergence. In Tab. [2] the exact resonances computed this way are 
compared to the corresponding semiclassical values. The relative error reaches 
its maximum value of about 5 percent at the first state |0, 0; +, +). The relative 
error shrinks rapidly for larger real parts of the resonance energies. 

As shown in Fig.[l3]the resonance energies (resp. the corresponding wavenum- 
bers) form a grid in the complex energy (resp. wavenumber) plane. The grid 
sites can be labeled by the quantum numbers and and the parities tt^ and 
TTy. As also shown in Fig. [HI each step of the cumulative reaction probability 
is associated with one 'string' of resonances of fixed and iTy. Wc note that 
interestingly the values of S2 are located along a smooth line in the complex 
S2 plane. This can be understood semiclassically from taking the quotient of 
the EBK quantized actions and I^; whose integrals both scale with Ve, and 
hence leads to the energy independent condition 

n4 + i(2-vr.)_ 4(^-^1)^ ^^^^^ 



on 82- 

In Fig.[TH]we present the separated wavefunctions ipi, and tpa- for a selection of 
resonances. We see that the scaled wavefunction obey the boundary conditions 
determined by the parities tTx and tTj, and the exponential decay as u — s- cx), and 
also have the expected number of nodes determined by the quantum numbers 
n,^ and n^. 

For a selection of resonances, the total density of position given by 

|V'(i.,A)|2 = |V'.WA(A)P (138) 

is shown in Fig. [11] Note that there are no nodal lines other than the coordinate 
axes for states with negative parities. This is due to the complex valuedness 
of the energy E and the separation constant for these states. In particular 
sections of these plots along the bottleneck x = (which would give the 2D 
analogue of Fig.[2T]in the 3D case) would lead to densities that are greater than 
zero at every y G (— 1, 1) apart from a possible zero at ?/ = if Tr.y = — . 



6.2 The 3D system 

6.2.1 Semiclassical computation of resonances 

In the 3D case wc combine the complex EBK quantization condition p28p which 
we write in terms of the symmetry reduced action analogously to (I130p with 
the uniform quantization conditions (|112p . The resulting set of equations 

ic^{nc + ^)h, i^ = {n^ + ^)h, = -i(n5 + i(2-7r,))?i, (139) 

where n(^,nr,,n(^ G INg, and and are defined in (|113p and (|114p . can again 
be solved by a Newton procedure on the (real) six-dimensional space of complex 
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Figure 19: Two-dimensional contourplots of the densities of position in the x — y 
plane for all the resonance wavefunctions p38p of the 2D system, with = 5, 
for all the resonance states shown in Fig. dH Light blue corresponds to low 
probability, red corresponds to high probability. 



Si, S2 and E. To this end the integration paths defining the integrals that enter 
the quantization conditions (|139p have to be continued into the complex plane as 
shown in Fig. fOT d.e). We label the resulting resonance modes by the Dirac kets 
jn,^, 71,,, n^;Try,TTz, tTx)- For fixed quantum numbers n^, n,, and parities tTj^, tt^ we 
find the resonances by starting near to, but not at, the corresponding "opening" 
tick on the real energy axis of Fig. 1 151 and by smoothly moving the parameter 
into the complex plane. We then go through the grid {n^, tTx) G INq x { — 1, 1}. 
As with the 2D case, special care has to be taken of resonances which have 
energies close to the imaginary axis. We give a list of semiclassically computed 
resonances in Tab. [3l 



6.2.2 Exact computation of resonances 

For the computation of the exact quantum resonances using the complex scaling 
method we again transform the A component of the separated wave equation 
(see dSll)) 

-£S'=l<«*«-^«'<-^'+"*- <"») 

We are now looking for a transformation of A and analogous to (jl32p and 
()133|) in the 2D case which yields a system of type 'kinetic-plus-potential' with 
the potential going to zero at infinity. In the 3D case this can be achieved by 
setting 

(t(A) ^ag^tn(A,g) =ag^ sn(A, g) 

cn(A, q) 

and 

^,(A) = -4±L = ./4 ^ , (142) 

Vda/dA V 1' ((^2 _^a2)(^/2^2 _^a2))i/4 
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N ^ " '/ " t 








Iqu, 


^.-< 






A.E 





+ + + 


3.16SS - i0.5S35 


2.4498 + i0.4141 


0.4408 + i0.0771 


2.9886 - 10.5942 


2.5052 + 10.4464 


0.5802 + 10.1081 


5.4 





+ + 


2.8348 — il.7219 


2.0997 + 11.1533 


0.3756 + 10.2147 


2.6279 - 11.7526 


2.0976 + 11.2224 


0.4814 + 10.2966 


4.8 





+ + 


7.7066 - iO.9376 


2.5633 + i0.2801 


1.0338 + i0.1196 


7.5633 - 10.9414 


2.5844 + 10.2860 


1.0528 + iO.1237 


1.8 





+ 


7.3621 - 12.7939 


2.4011 + iO-8131 


0.9646 + iO-3473 


7.2031 - 12.8049 


2.4139 + iO.8284 


0.9790 + iO.3584 


1.7 





+ - + 


8.6233 - iO.9976 


2.6146 + iO-2627 


0.6062 + i0-0566 


8.4030 - i0.9971 


2.5386 + iO.2689 


0.6767 + i0.0661 


2.6 





+ 


8.2655 - 12.9751 


2.3693 + iO-7649 


0.4739 + iO-1647 


8.0398 - 12.9732 


2.3853 + i0.7816 


0.6391 + i0.1921 


2.4 


10 


+ + + 


14.2036 - il.2881 


2.6094 + i0-2104 


1.3141 + i0-1140 


14.0702 - il.2859 


2.6202 + iO.2116 


1.3142 + i0.1147 


0.93 


10 


+ + 


13.8402 - i3.8499 


2.6176 + iO-6197 


1.2643 + iO-3359 


13.7074 - 13.8435 


2.5266 + i0.6231 


1.2634 + iO.3376 


0.90 





+ 


14.7294 - 11.3158 


2.5926 + i0-2060 


1.1386 + i0-0981 


14.6744 - il.3161 


2.6009 + 10.2065 


1.1176 + i0.0063 


0.37 







14.3639 - 13.9332 


2.6040 + i0.6072 


1.0964 + iO-2892 


14.3106 - 13.0344 


2.5112 + 10.6083 


1.0767 + iO.2837 


0.34 


10 


+ + + 


17.6729 - 11.4364 


2.6143 + i0.1856 


0.3219 + i0.0250 


17.3973 - il.4300 


2.5255 + 10.1866 


0.3476 + i0.02 70 


1.0 


10 


+ + 


17.2100 — 14.2966 


2.4407 + iO.5486 


0.3120 + i0.0738 


17.0367 - 14.2775 


2.4507 + 10.5510 


0.3368 + 10.0799 


0.97 




















7 3 


+ + + 


514.6176 - 17.8748 


2.6479 + i0-0357 


1.4866 + i0.0219 


514.4407 - 17.7908 


2.6483 + 10.0363 


1.4868 + i0.0217 


0.04 


7 3 


+ + 


514.2416 - 123.6217 


2.6463 + i0.1069 


1.4840 + i0.0656 


514.0764 - 123.3700 


2.6457 + 10.1068 


1.4842 + i0.0650 


0.03 


6 3 


+ 


514.6179 - 17.8748 


2.6479 + i0.0357 


1.4856 + i0.0219 


514.4407 - 17.7908 


2.6483 + 10.0353 


1.4858 + 10.0217 


0.04 


6 3 




514.2419 — 123.6217 


2.6453 + 10.1069 


1.4840 + 10.0656 


514.0764 - 123.3701 


2.6457 + 10.1058 


1.4842 + 10.0650 


0.04 


8 2 


+ + 


514.8034 - 17.9062 


2.6919 + 10.0360 


1.9256 + 10.0287 


514.5142 - 17.8196 


2.6923 + 10.0356 


1.9263 + 10.0284 


0.06 


8 2 


+ 


514.4227 — 123.7160 


2.6893 + 10.1080 


1.9234 + 10.0859 


514.1450 — 123.4562 


2.6897 + 10.1069 


1.9242 + 10.0851 


0.06 


8 2 


+ - + 


514.8034 - 17.9062 


2.6919 + 10.0360 


1.9256 + 10.0287 


514.5142 - 17.8196 


2.6923 + 10.0356 


1.9263 + 10.0284 


0.06 


8 2 


+ 


514.4227 — 123.7160 


2.6893 + 10.1080 


1.9234 + 10.0859 


514.1450 — 123.4562 


2.6897 + 10.1069 


1.9242 + 10.0851 


0.06 


2 7 


+ + + 


515.1296 - 17.8539 


2.5566 + 10.0349 


0.5719 + i0.0083 


514.9221 — 17.7729 


2.5570 + 10.0346 


0.5730 + 10.0082 


0.04 


2 7 


+ + 


514.7588 - 123.5592 


2.5540 + 10.1047 


0.5713 + i0.0249 


514.5623 — 123.3165 


2.5544 + i0.1036 


0.6723 + 10.0246 


0.04 


4 6 


+ + 


518.0938 - 17.8509 


2.5998 + i0.0351 


1.0038 + i0.0143 


518.0400 - 17.7650 


2.6000 + i0.0348 


1.0033 + i0.0141 


0.01 


4 5 


+ 


517.7246 - 123.6602 


2.6972 + i0.1053 


1.0028 + i0.0429 


517.6817 - 123.2929 


2.5974 + i0.1042 


1.0022 + i0.0423 


0.03 


8 


+ - + 


519.6687 - 17.8903 


2.6062 + iO-0344 


0.0680 + iO.OOlO 


519.4248 - i7.8110 


2.5066 + i0.0340 


0.0688 + iO.OOlO 


0.01 


8 


+ 


519.1991 - 123.6686 


2.6036 + i0.1031 


0.0680 + 10.0030 


519.0642 - 123.4308 


2.5040 + 10.1021 


0.0687 + 10.0030 


0.03 



Table 3: The quantum mechanical complex eigenvalues {Eq^n^kq^^^qm) and the semiclassical complex eigenvalues {iJscj ^sc, ^sc} of the 
3D asymmetric hyperboloidal billiard for the ranges Rc(i?) < 20 and 500 < Re(£') < 520 and for {nj,7rx} = {0,±} (the first two families 
of resonances). The relative error t^E = {\Eqm\ — \Esc\)/\Eqm\ is given in percent, (h = 1, m = 1.) 



where the modulus of the eUiptic functions is again given hy q = c/a and 
q' = {1 — (7^)^/^. This leads to the new wave equation 

" £^ = (i^.,eff - K,cff)^. , (143) 

with the effective energy and potential given by 

E[{2k - a^q^){q'^a^ + a^) - l]/q'^ + h\6q'^cr^ + + q'^)/4m 



[a"^ + a'^){q''^a'^ + a^) 
3(2(7'V3 + (a2+g'2)^)2 



2m4(CT2 + a2)2(q/2^2 _^a2)2 



(144) 



respectively. 

As opposed to the effective potential in l|140p . the potential T4,cff goes to 
zero as \a\ 00, and accordingly, for real E, and S2, the wave function ipa- is 
a plane wave as u ^ 00. In fact, in terms of the original Cartesian coordinates 
{x,y,z) the solutions have to become plane waves for 00. From ((20|) and 

noting that ^ = adn(A, g)/ cn{X,q) we see that a = aq'tii{X,q) is proportional 
to X and this is the motivation for the transformation p4ip . The scaling of 
the wavefunction (|142p is performed in order to again obtain a system of type 
'kinetic-plus-potential' . 

Similarly to the 2D case, to compute the resonances we substitute a in (jl43p 
with ere'". We have to solve the C and rj components of the wave equation and 
the corresponding boundary conditions described in Sec. 15.2.11 in combination 
with the equation for ip^ in (|143p with the boundary conditions at zero deter- 
mined by the parity tTx and V'ct(o') — > as |cr| — > cx) (see the analogous conditions 
for the 2D case in Sec. 16.1."^ . In our numerical procedure, which consists of a 
shooting method like in the 2D case, we again choose a sufficiently large CToo at 
which we require '0<t(o'oo) = 0. Decomposing all equations with respect to their 
real and imaginary parts the shooting method results in a (real) six-dimensional 
Newton procedure acting on the complex three-dimensional plane E — S1 — S2 (or 
equivalently E—k — l). Using the semiclassical values from the previous section 
as the starting values the shooting method always converges to the expected res- 
onance state. Like in the 2D case, special care has to be taken of those resonance 
states of energy close the imaginary axis. In Tab. [3] the exact resonances are 
compared to the semiclassically computed resonances. The relative error of the 
semiclassical complex energy eigenvalues reaches its maximum value of about 
5 percent for the first state |0, 0, 0; -|-, -|~, The relative error shrinks rapidly 
for larger resonance energies. Plotting the resonances in the complex energy 
(resp. wavenumber) plane in Fig. [15] we see that each string of resonances of 
fixed quantum numbers tiq and n^, and parities tt^ and tt^ give rise to one step 
of the cumulative reaction probability. Since we have three quantum numbers 
for the 3D system (as opposed to the two quantum numbers in the 2D system) 
the resonances can be viewed to form the superposition of an infinite number 
of grids of the more regular type of grid found in the 2D system in Fig. [T31 

From studying the separated transverse wavefunctions ipi, and ip^j, one can 
see that their real and imaginary parts are similar, but not equal. Also for 
increasing quantum number n\ these wavefunctions vary little (see Fig. 1201 ) 



52 





1.0 20 



Figure 20: Real parts of the separated resonance wavefunctions (solid lines) 
of the eight states |n^, n^, n^; %, tTz, t^x) — |0, 2, {0, 1}; +, — , ±) (classically BB, 
upper panels) and \n(^, rirj, n^; %,7rz, tt^;) = |4, 0, {0, 1}; +, +, ±) (classically WG, 
lower panels). For fixed {ri|j, n^, tTj,, tt^} (and hence fixed {rij^jn^}), ijj^, and i/i^ 
are the same (far left panels) whilst for increasing n\ the center panels show 
the increase in nodes of the separated wavefunction in the reaction coordinate 
(J. The dashed lines are the real parts of the corresponding effective potentials 
Vs_csis)- The wavefunctions are plotted at the real part of the effective energies 
Eg^eS (dotted lines). The two far right panels are the real parts of two sample 
effective potentials V^,eff viewed at long range. 
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Figure 21: Probability contours in the section a; = of the wave functions 
of the resonance modes |n^, n^, n^; tTj,, tt^, tt^;). These correspond to the first 
three famihes of resonances associated with the steps labelled |7) to |13) in 
Fig. [jni For TTx ^ — the ellipses shown here are actually nodal surfaces of 
the total wavefunction. The separated wavefunction is therefore ignored 
for the purposes of this figure. Light blue corresponds to low probability, red 
corresponds to high probability. Only axis lines (i.e. the symmetry lines) are 
nodal. 

However, the real and imaginary parts of -tpi, and -0^ have their zeros at slightly 
different values of the respective coordinates v and fi. Due to the complex 
valuedness of and s| the total wavefunctions do not have any nodal surfaces 
apart from the ones along the symmetry planes if the corresponding parity is 
negative. This can also be seen in Fig. [^11 where we present the contours of the 
resonance wavefunctions as their intersections with the y — z plane (ignoring 
the component ipx of the total wave function which would lead to zero valued 
total wavefunctions when tt^, — —). Note that these contours look quite different 
from the analogous contours for the scattering states in Fig. [T71 In addition to 
the absence of nodal lines there is also (as to be expected) no clear distinction 
between whispering gallery and bouncing ball modes. 
The total probability density is given by 

\^{y, ^l, A)p - \Ml^)M^^)M^)? ■ (i45) 

In order to compute the separated wavefunction ■(/'a (A) we use again the trans- 
formation p42p and the inverse of (|14ip which is given by the elliptic integral 

/•cr 

X(a)= tn-\<7/aq',q)= ^ ^^^^ ■ (^46) 

In Fig. [21] we show the isosurfaces of (|145p inside the hyperboloidal boundary for 
two examples of resonance states which further illustrates the absence of nodal 
surfaces (in addition to the Cartesian coordinate planes for negative parities). 

7 Conclusions and Outlook 

In this paper we demonstrated how ideas from transition state theory can be 
used to compute the classical and quantum mechanical transmission proba- 
bilities for transport through entropic barriers. For barriers associated with 
saddle points of the potential (and more generally saddle type equilibria of the 
Hamilton function) it has recently be shown that the transport through the 
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Figure 22: Three-dimensional contourplots of the probabihy densities of the 
wavefunctions of the 3D system, with parameters (? — 0.2 and a} — 5, 
for the two resonance states n^, rt^; tTj,, tt^, tTj;) = |0, 0, 0; +, +, +) (left) 
and jn^;, n^, n^; TTj^, TTx, TT^r) ~ |1, 1, 0; +, +, — ) (right). The isosurfaces are 
|V(J^,M,A)|2 = 0.1 (cyan), \'>\^(y,p.,\)\^ = 0.3 (green), |V'(i^, A)|2 = 0.5 (red) 
and |V'('^, A)P = 0.9 (blue). Note that there are no nodal surfaces apart from 
the y — z plane due to a negative parity tt^ for the second state. 
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phase space bottlenecks induced by such saddle points are controlled by a set 

of phase space structures. In the present paper we identified the analogous 
phase space structures for potentialless barriers where the phase space bottle- 
necks are induced by hard wall constrictions. We focused on the special case of 
hyperboloidal constrictions in two and three dimensions for which the classical 
and quantum transmission problems are separable, and hence facilitate a very 
detailed, and to a large extent also analytical, study. For these systems, we 
showed that like in the case of smooth systems one can construct a dividing 
surface which has the property that it is crossed exactly once by reactive tra- 
jectories, and not crossed at all by nonreactive trajectories. This, like in the 
smooth case, leads to a rigorous realization of Wigner's transition state theory, 
i.e. to an exact computation of the classical transmission probability from the 
(directional) flux through the dividing surface. Similarly the quantum cumula- 
tive reaction probability can be computed exactly from the quantum mechanical 
flux through the dividing surface. We showed that like in the case of smooth 
systems the dividing surface is linked to an unstable invariant subsystem with 
one degree of freedom less than the full system. In the context of chemical 
reactions such an invariant subsystem located between reactants and products 
is referred to as the transition state or the activated complex. For the 2D and 
3D systems studied in this paper the transition states consist of the billiard in a 
one-dimensional square well and the billiard in an ellipse, respectively. Like in 
the smooth case the transition states in the hard wall constrictions have stable 
and unstable manifolds which have sufficient dimensionality to form separatrices 
which separate reactive trajectories from nonreactive trajectories and thus play 
a key role for the classical transmission problems. 

Quantum mechanically the transition states manifest themselves on the one 
hand through a 'quantization' of the cumulative reaction probability and quan- 
tum resonances. The quantization of the cumulative reaction probability refers 
to the stepwise increase of the cumulative reaction probability as a function 
of energy each time a new state fits into the respective transition state. In 
fact, as discussed in some detail, the cumulative reaction probability is approx- 
imately given by the integrated density of states of the invariant subsystems 
associated with the transition states. For the quantum systems, the Hcisenberg 
uncertainty relation excludes the existence of invariant subsystems analogous to 
the classical case. Instead a wavcpackct initialised on the classically invariant 
subsystems will decay (exponentially fast) in time with the lifetimes being de- 
scribed by the resonances. We computed such resonances semiclassically from 
the poles of a meromorphic continuation of the semiclassical expression for the 
transmission probability to the lower half of the complex energy plane. We 
showed that this leads to a very good agreement with the results obtained from 
the numerical computation of the exact resonances from the complex scaling 
method. The separability of the systems yields an assignment of the resonances 
by quantum numbers. We showed that each string of resonances corresponding 
to fixed quantum numbers associated with the transverse degrees of freedoms 
gives rise to a step of unit size of the cumulative reaction probability. Despite of 
their separability, the systc;nis studied in this paper display quite a rich variety 
of dynamics. In particular the transition state of the 3D system involves two 
different types of modes that we referred to as whispering gallery and bouncing 
ball modes. We showed that the energetic quasidegeneracy of the whispering 
gallery modes leads to steps in the cumulative reaction probability of effective 
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stepsize 2. Similarly the corresponding resonances are quasidegenerate with 
respect to their complex energies. 

For many aspects of our study we heavily used the separability of the sys- 
tems. We in particular presented a detailed study of the analytic nature of the 
integrals which enter the EBK and uniform semiclassical quantization schemes 
as action and tunnel integrals. We showed that all these integrals can be in- 
terpreted as Abelian integrals on an elliptic (2D system) or hyperelliptic (3D 
system) curve with branch points determined by the separation constants and 
the geometry parameters of the hyperboloidal constrictions. This interpretation 
led to a natural extension of the semiclassical quantization conditions of scat- 
tering states to quantization conditions for resonances states by complexifying 
the separation constants and integration paths. 

Though the systems discussed in this paper are special due their separability 
many of the phase space structures found are expected to exist in systems with 
more general constrictions that lead to nonseparable dynamics. In particular, 
it is to be expected that the phase space structures associated with the transi- 
tion states persist under (small) deformations of the hyperboloidal constrictions 
which (generically) would destroy the integrability of the systems. For smooth 
systems, this has been studied already in quite great detail. The phase space 
structures for such systems can, e.g., be practically determined from a normal 
form expansion about the saddle equilibrium points [U [8] . For systems where 
the barriers are induced by hard wall constrictions the determination of the 
phase space structures associated with the transition state is more challeng- 
ing in the generic nonseparable case. The best approach here seems to be the 
study of these structures in terms of a billiard map, i.e. a map resulting from 
taking snapshots of the billiard dynamics from one specular reflection to the 
next. For a 2D system the resulting map is a symplectic map from to E,^, 
with the periodic orbit that forms the transition state in this case appearing 
as a hyperbolic fixed point of the map. There are well established methods for 
determining the fixed point and also (using methods based again on a normal 
form) for computing the stable and unstable manifolds of the fixed point. How- 
ever, the situation is much more involved in the 3D case, where the resulting 
billiard map is a symplectic map from to R^. The transition state then forms 
a two-dimensional manifold in the domain and image of this map. Its stable 
and unstable manifolds arc three-dimensional. Computing these manifolds in 
practice is quite difficult. This gives an interesting field for future studies. 
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